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We derive the gravitational radiation-reaction force modifying the Effective One Body (EOB) 
description of the conservative dynamics of binary systems. Our result is applicable to general 
orbits (elliptic or hyperbolic) and keeps terms of fractional second post-Newtonian order (but does 
not include tail effects). Our derivation of radiation-reaction is based on a new way of requiring 
energy and angular momentum balance. We give several applications of our results, notably the 
value of the (minimal) "Schott" contribution to the energy, the radial component of the radiation- 
reaction force, and the radiative contribution to the angle of scattering during hyperbolic encounters. 
We present also new results about the conservative relativistic dynamics of hyperbolic motions. 
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I. INTRODUCTION 

The Effective One Body (EOB) formalism [H-H is 
an approach to the relativistic dynamics of gravitation- 
ally interacting binary systems which was originally pro- 
posed as a way to extend the validity of the usual 
post-Newtonian (PN) formalism beyond the slow-motion 
(v 2 /c 2 < 1) and weak-field (GM/(c 2 r) < 1) regime. The 
EOB approach is made of three, basic building blocks: 

1. a description of the conservative (Hamiltonian) 
part of the dynamics of two compact bodies; 

2. an expression for the radiation-reaction force T 
which must be added to the conservative, Hamil- 
tonian equations of motion; 

3. a description of the asymptotic gravitational wave- 
form emitted by the binary system. 

The building block 1, i.e., the EOB Hamiltonian, has 
been analytically computed with an increasing accuracy 
in a sequence of papers, both for non-spinning black holes 
0, Q , for spinning black holes 048 1 and for systems in- 
volving tidally-deformed bodies 0, [l(| ■ In addition, the 
comparison between the EOB dynamics and numerical 
simulations of binary systems has allowed one to improve 
the knowledge of some of the functions entering the EOB 
Hamiltonian (see Ref. [ll[ for a review). More recently, 
results from gravitational self-force theory [l2| have also 
allowed one to learn new information about the EOB 
formalism (See Ref. [l3| for recent progress and refer- 
ences). The description of the second building block, i.e. 
the radiation-reaction force T has also improved over the 
years, both through the conception of new resummation 
methods [14| and from the comparison with numerical 
simulations (both in the comparable- mass case [HI, [l6[ , 
and in the extreme-mass-ratio case [1343 )■ The same 



remarks apply to the third building block, i.e., the grav- 
itational waveform. 

While the EOB Hamiltonian is able to describe the 
conservative dynamics of general binary orbits (quasi- 
circular, elliptic-like or hyperbolic-like), the currently ex- 
isting accurate implementations of the radiation-reaction 
force and of the emitted waveform are limited to the case 
of quasi-circular, inspiralling orbits. The main reason be- 
hind this limitation is that the EOB program was orig- 
inally motivated as a tool for computing accurate wave- 
forms from the type of circularized binary systems that 
are likely sources for ground-based interferometric gravi- 
tational wave detectors. However, the progress in numer- 
ical relativity simulations has opened the possibility of 
numerically exploring the dynamics of binary systems in 
more exotic configurations. For instance, Refs. 0, l2l| 
have considered high-velocity encounters of black holes 
and other bodies, and Ref. [22| has considered eccentric 
orbits of binary black holes. We anticipate that more 
simulations of general orbits will become routinely possi- 
ble in the near future. See Ref. [23| for a recent example, 
and more references. 

This perspective motivates the main aim of the present 
work, namely, to provide an expression of the radiation- 
reaction force jF along general orbits (elliptic or hyper- 
bolic) within the EOB formalism. [We leave to future 
work a corresponding generalization of the EOB gravita- 
tional waveform.] 

Gravitational radiation-reaction, notably in binary 
systems, has a long history. Let us only recall that three 
general different approaches have been used. The first 
approach derives the full equations of motion of matter 
(including both conservative and radiative effects) from 
a direct integration of the retarded field generated by the 
source. Because of its difficulty, this approach has been 
implemented essentially only up to the next-to-leading 
order in T , i.e., at the fractional 1PN accuracy [23 - |28j . 
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A second approach focuses on the radiation-reaction 
piece in the equations of motion and derives it by us- 
ing a matching between between the near-zone field and 
the wave-zone field. This approach has been also imple- 
mented only up to the next-to- leading order in T |29l - l38j , 
with some vistas on the effect of tails (39J. 

Finally, a third approach is based on requiring a bal- 
ance between the losses of mechanical energy and an- 
gular momentum radiated by gravitational waves at in- 
finity. This "third" balance approach has been particu- 
larly developed by Iyer and Will and their collaborators 
[IS H3i E3| and has been implemented to a higher PN 
accuracy than the other approaches, namely the next- 
to- next-to- leading order in IF, i.e., the fractional 2PN 
accuracy [40j |. 

Note, however, that Ref. [40] does not include the 
effect of tails. We shall, similarly, postpone the inclusion 
of tails (entering at the fractional 1.5PN, w 3 /c 3 , level) to 
future work. We note that Ref. [3!| has shown that the 
tail contribution to IF satisfies the balance requirement. 

In view of the technical efficiency of the balance ap- 
proach (and of the direct proof by several authors of 
the consistency between this approach and other ones 
[33l HH), we shall also base our work on this ap- 
proach. However, instead of attempting to "translate" 
the radiation-reaction force T derived in Refs. [33l . l34l . l40j 
(which was derived in harmonic coordinates, and was ex- 
pressed in terms of quasi-Newtonian equations of mo- 
tion) into the EOB formalism (which uses different coor- 
dinates, and Hamiltonian equations of motion), we found 
more efficient to develop a new way of using the balance 
approach. We shall explain in detail below our new way 
of implementing the balance approach. 

Let us only say here that it is based on three essential 
ingredients: (i) we start from the 2PN-accurate expres- 
sions of the fluxes of energy and angular momentum, $ e 
and 4>j, that have been derived in the PN literature (4lT - 
l45j (see references [46M50lj for recent higher PN accuracy 
results). These fluxes are expressed in terms of three 
scalars v 2 ,, f\ and GM/r^, where X/j and v/,. denote har- 
monic coordinate and velocities (of the relative orbit). 
Then, (ii) we derive the transformation connecting the 
three scalars v^, r\ and GM/r^ to the three scalars that 
are natural within the EOB formalism, namely p e 2 , p\ r 
and GM/r e , where x e and p e denote EOB coordinates 
and momenta. Finally, (iii) we introduce a new way 
of using the two EOB-expressed fluxes $s(x e ,p e ) and 
<I>j(x e ,p e ) to derive the two independent components of 



the radiation-reaction force Fr B (x e , p e ) = IF' 
J-( cob) (x e ,p e ) = (x e x^ cob )).e z . 

The structure of this paper is as follows. We present 
in Sec. II our new way of implementing the balance ap- 
proach. Then, in Sec. Ill, after presenting a brief review 
of the EOB formalism, we apply our method to the 2PN- 
accurate EOB-variables forms of &e and $j, and derive 
explicit expressions for IFr" ^ and J-^° h \ We also ob- 
tain the explicit expressions of the associated "Schott" 



_ T-(eob) 



n, 



energy contribution. Sec. IV discusses the gauge free- 
dom in IF and explains how it is related to the freedom 
in defining the Schott contributions to the energy and 
angular momentum. Then, Sec. V gives some applica- 
tions of our results, and discusses notably the scattering 
angle during hyperbolic encounters, and its modification 
by radiation-reaction effects. We summarize our main 
results in Sec. VI, and discuss future directions. Finally, 
to relieve the tedium we have relegated several explicit 
technical details to various appendices. 



II. A NEW APPROACH TO 
RADIATION-REACTION 

Here, we introduce a new approach to the computation 
of radiation reaction by the balance method. Let us con- 
sider the effect of adding a radiation-reaction force, say 
IFi, to the Hamiltonian form of the equations describing 
the relative motion of a binary system (with masses mi 
and 7712) 



d^(x,p) 

dpi 



Pi 



dx l 



IFi 



(2.1) 



Here "H(x, p) denotes the Hamiltonian and a dot denotes 
differentiation with respect to time. When considering 
the motion within the orbital plane, we can take as co- 
ordinate and momenta x l = (r,<p) and pi = (p r ,p<j,). 
Correspondingly, the radiation-reaction will have two in- 
dependent components: J> and 

Let us see how one can determine the two force com- 
ponents T r and F$ by writing balance equations for the 
energy and the angular momentum of the binary system, 
namely 

^(system) (*) = H(x(t),p(t)) - (mi + m 2 )c 2 

■/(system) (*) = P<f>(t) ■ (2-2) 

On the one hand, the equations of motion (|2.1|) 
yield the following time changes for E( system )(t) and 



J, 



(system) 



(0 



dE, 



(systc 



dt 



)(t) _du _ ..m . &r_ 

j, x n : T Pi n X J~i 

dt ox 1 dpi 



d J [system) ft) _ dp^ _ 8% 

dt dt d<j) 



Frh 



(2.3) 



The explicit form of these two equations read (when using 
the fact that T-L does not depend on tp) 



-^(system) (*) — r ^r + fy? <j> 



and 



(system) 



it) 



~dt 



— F,h 



(2.4) 



(2.5) 



It will also be useful to consider the following combina- 
tion of these two equations 



-^(system) (system) rF r . 



(2.6) 
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Formally speaking, Eqs. (|2.4[) and (|2.5[) provide two 
equations relating the two unknowns T T and T$ to the 
losses of energy and angular momentum. 

On the other hand, we require that there is a balance 
between the energy and angular momentum losses of the 
system, and the corresponding energy and angular mo- 
mentum fluxes (in the form of gravitational waves) at in- 
finity, say $e and $,/. As was pointed out by Schott long 
ago [51|, one cannot, however, simply equate -EVsystem) 
and J( S y S tem) to, respectively, — and — $j. One must 
allow for the existence of Schott terms that represent 
additional contributions to the energy and angular mo- 
mentum of the system, due to its interaction with the 
radiation field, say 2£( sc hott) (*) = ^(schott) p(*)) and 
■/(schott) (*) = -/(schott) (*(*),p(t)). The correspondingly 
modified balance equations then read 



E, 



(system) 



E, 



(schott) 







'(system) + -/(schott) + $J = . (2.7) 

Inserting the identities ([13}, (j2"?5|) into (pTT)) leads to the 
following two conditions on the two components of the 
radiation-reaction force 



•F 6 + J( 



E, 



(schott) 



(schott) 







J 



0. 



(2.8) 



Up to now, all the equations we have written down 
are equivalent to the standard "balance approach" to 
radiation-reaction, as used, in p articular, by Iyer, Will 
and collaborators [33|, H3, So, Ef , except for the fact that 
we are working within a Hamiltonian framework. Let 
us now explain the new, simplifying features of our ap- 
proach. 

The first simplifying feature is to note that it is al- 
ways possible to impose the condition that the Schott- 
contribution to the angular momentum vanishes: 



-/(schott) (x(f),p(f)) = 0. 



(2.9) 



The proof that this is possible is simply that, after im- 
posing Eq. (|2.9p . we shall be able to find a solution to 
the general balance equations (|2.8|) . Indeed, after making 
the assumption (|2.9j) . we can use the second Eq. (|2.8|) to 
determine the instantaneous value of the (^-component of 
the radiation-reaction force, in terms of the correspond- 
ing instantaneous J-flux: 



^ = -$j(x(t),p(t)). 



(2.10) 



Let us note in passing that the result (|2.10[) for T$ 
is standardly used in the current implementations of the 
EOB equations of motion [TTJ] . Then, by inserting the 
result (|2.10[) into the first equation (|2.8|) , we get an equa- 
tion involving only JF r and i?( sc hott) > namely 



rT r + E. 



(schott) 



= -$ 



EJ > 



(2.11) 



where we introduced the notation 

$ EJ (x ) p) = $ E (x ) p)-0(x,p)$ / (x,p). (2.12) 



As we shall discuss in detail in the next section, we as- 
sume here that we have in hands explicit expressions 
for $b, $,/ (as well as for the "combined flux" §ej) 
as functions of the instantaneous dynamical state of 
the system. Within a Hamiltonian framework it means 
$ B = $ B (x,p), $j = $,/(x,p) and $> E j = $ B j(x,p). 
[Note that, by Hamilton's equations, the instantaneous 
orbital frequency 4> entering $ e.j is a function of position 
and momenta, given by </>(x, p) = <9"H(x, p)/dp c f ) . As we 
shall further discuss below, contrary to $b and $j, is 
not a gauge invariant quantity; we shall only consider its 
explicit expression 0(x, p) in EOB coordinates.] 

While Eq. (|2.10[) provides an explicit expression for 
in terms of the instantaneous state of the system, our 
remaining problem is to show how Eq. (|2.11l) can be 
used to determine both J>(x, p) and /^(schott) ( x , p)- Let 
us now explain how this can be done. 

The basic idea is that the specific combination $ e j has 
the property of vanishing along circular motions. Indeed, 
it is well known that (because of the monochromatic na- 
ture of the emitted radiation) one has = il&j along 
a circular motion with orbital frequency f2. As a conse- 
quence, when considering general, noncircular motions, 
$B7 will necessarily be given by an expression which can 
be written as a combination of the two independent quan- 
tities that vanish along circular motions, namely 



Zi(x,p) =vl 



and 



Z 2 (x,p) 



dr 



-rp r 



(2.13) 



(2.14) 



where the factor r in Z 2 is introduced for later conve- 
nience. [See Sec. IIIC where we will work with rescaled 
versions of Z\ and Z-i that have the same dimensions.] 

Here, we are availing ourselves of several simplifica- 
tions that are allowed at the PN accuracy at which we 
shall be working. First, as we shall explicitly check, the 
combination $£j(x, p) is invariant under time reversal, 
and can therefore be expressed as a function of <~ r 2 , 
rather than simply of p r ~ r. Second, modulo terms of 
5PN order (i.e., 0(l/c 10 )) one can neglect the J> contri- 
bution to the link between p r and —8% /dr. 

We can then write 



$ £J (x,p) = $i(x,p)Zi(x,p) + $ 2 (x,p)Z 2 ( x ,p) 

(2.15) 



-r 2 ^ dp r 



where <&i and $2 exist but are not uniquely defined. For 
instance, we can move a term oc Z2 in $1 to $2, and 
reciprocally a term cx Z\ in <I>2 to $1. We shall discuss 
below the effect of these ambiguities in the definition of 
$1 and <j> 2 . 

Operating by parts on the second expression (|2 . 1 5[) 
(which involves p r ), we can then write 



$£j(x,p) =Pr 



Pr<S>l 
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(2.16) 



which is a decomposition of <&ej in a part proportional 
to p r (and therefore to r, in view of r = dH/dp r ), and a 
total derivative. But, such a decomposition is precisely 
the content of the balance requirement (|2.11l) . 

We therefore see that, given any choice of $1 and $2 
such that Eq. (|2.15[) holds, we can obtain one particular 
corresponding solution to Eq. (|2. 1 1|) . namely 



JV(x,p) 



Pr_ 

r 



dt 



(r$ 2 ) 



-B(schott) (x, p) = rp r ® 2 



(2.17) 



In keeping with our approximations, the time deriva- 
tive of r$2 (x, p) in the first Eq. (|2.17j) should be evalu- 
ated along the (conservative) Hamiltonian dynamics, so 
that T r can be explicitly expressed in terms of the in- 
stantaneous dynamical state of the system. 

The results ([2~T7| . together with Eqs. fl21JJ and (j2TTUl) . 
give a constructive algorithm for determining the two 
components T r and of radiation-reaction, as well 
as the Schott contributions to energy and angular mo- 
mentum. [This contrasts with Refs. [33l . [34], |4^] which 
had to use the method of undetermined coefficients.] 
This proves our claim that is indeed possible to define 
a radiation-reaction force such that the Schott contribu- 
tion to the angular momentum vanishes. [By contrast, 
one can show that it is generally impossible to define 
T% such that -E( sc hott) vanishes.] We shall discuss later, 
while implementing our construction, the impact of the 
non uniqueness in the decomposition (|2.15[) , as well as a 
simple, algorithmic way of fixing it. Let us only note here 
that, in keeping with the analysis of Iyer and Will [33ll33| 
and later developments by Gopakumar et al [4(J, all the 
non uniqueness in the definition of the radiation-reaction 
F has the character of a gauge freedom (and is actually 
related to possible coordinate changes) . This also applies 
to the freedom of setting J( sc hott) to zero, that we have 
used here to simplify the search for T. 



III. RADIATION REACTION FORCE IN THE 
EOB FORMALISM 

Let us now apply the method explained in the previous 
section to the construction of the radiation-reaction force 
in the EOB formalism. To do that, we need the following 
items: 

1. The expressions of the various flux functions $e, 
<&j and <&ej in terms of the positions and momenta 
of the EOB formalism; 

2. An algorithmic way of decomposing $bj(x, p) in 
the form (j2~T5]l . 

Before considering these items, let us recall the structure 
we shall need of the EOB formalism. 



A. EOB formalism: a short review 

At the 2PN accuracy that we shall consider here, the 
EOB Hamiltonian for the relative dynamics of two masses 
mi and rri2, is completely described by the following ef- 
fective metric 



ds 2 



= -A{r)c 2 dt 2 + B(r)dr 2 
+r 2 {d9 2 + sin 2 6»# 2 ) 



(3.1) 



where 



A(r) = 1-2 

B(r) = 1 + 2 

A(r)B{r) = D(r) 
Our notation is 



/ GM\ ( GM \ 3 

\ c 2 r J \ c 2 r J 

( GM \ 



\ c 2 r ) 



2(2-3 



\ c 2 r J 



+ 



, „ / GM\ 2 



n,r , m im2 fl 
M = Wil + 777,9 , = , V = 

It will often be convenient to work with 

GM 
u = . 



(3.2) 



(3.3) 



(3.4) 



With an abuse of notation we will then write 

u u 
A(u) = l- 2? +2^ + ... 

9 

11 11 

B{u) = 1 + 2— + 2(2 - 3^)— + ... 



D(u) = 1-6i/-j 



(3.5) 



The EOB Hamiltonian 'H( co b) is then defined as the fol- 
lowing function of the EOB coordinates (r, <j)) and mo- 
menta (pnPip) in the plane of the relative trajectory 



H 



(eob) 



Mc 2 J 1 + 2/. 



H 



(eff) 



where 



n 



(eff) 



= Mc 2 h . 

-.A(r) 
-A(r) 



(3.6) 



(n e ■ p e ) 2 , (n e x p e 
y?c 2 B{r) 



fj?c 2 



r,2 



c 2 B(r) c 2 r 2 



(3.7) 



that is 



(off) 



/-i<r 



A(u) 1 + 



A{u) 1 + 



A{u)p 2 r 


Pi 


D{u)c 2 


c 2 r 2 


A{u)p 2 , 


2 -2 

v> 3 


D{u)c 2 


c 2 



(3.8) 



5 



and 




(3.9) 



Here we have introduced a tilde to denote the result of a 
rescaling by the reduced mass fi, e.g. p = pj fj, and 



E : 



(eob) 



Mc 2 



(3.10) 



where we subtract the rest mass contribution to the en- 
ergy before scaling by fi. In addition it is convenient to 
introduce a special notation for some useful rescalings by 
GM, namely 



J = 



_ P<p 



GM fxGM 



GM 



t 



t 



GM 



(3.11) 



If we denote by V any quantity having the dimension 
of a velocity, we note that the dimensions of the GM- 
rescaled quantities u, j, q and t is u ~ V 2 , j ~ V^ 1 , 
q = |q| = it -1 ~ V~ 2 and i ~ V~ 3 . In the following, we 
shall often find convenient to work with the Hamiltonian 
pair of variables q,p r ] 4>,j- These variables are canon- 
ically conjugated with respect to the /i-scaled Hamilto- 
nian %(eob) = %(cot>)/Mj and correspond to an evolution 
with respect to the GM-scaled time t. For instance, we 
have 



GMT r 



dq _ dH(cob) dp r 5K(eob) 

di dp r ' di dq 

Ti =GM A=^T' s =J *- (3 - 12) 



Note also the vectorial relation 

r p J 

j = qx p = x - = 



GM [i GMfx 



(3.13) 



where J = r x p is the orbital angular momentum of the 
system. Let us also note the following relations 



an 



(cob) 



with 



dp a 



dpr 



1 cW( cff ) 
h Dpa 



c 2 A 



BH 



-Pr 



(off) 



r 2 H 



(3.15) 



(off) 



B. $ B , $j and <f> EJ in EOB variables 

Let us now indicate how one can express the flux 
functions &e, &j and $ej in terms of EOB variables. 



The first, crucial remark is that and $j are gauge- 
invariant quantities, and are scalars. [Note, however, 
that this is not true for <&e.j = &e — <f>$j, because 
<f> is not a gauge invariant quantity (along non-circular 
orbits), but depends on the chosen coordinate system. 
Here, we shall only consider the value of the combined 
flux in EOB coordinates: $|^ B = $ B - <p EOB <5>j] This 
implies that the numerical values of $ B and <!>j are inde- 
pendent of the choice of coordinates, and of any related 
choice of dynamical variables. We can therefore start 
from the results in the literature that have computed $ e 
and <&j, say at 2PN accuracy, in terms of, e.g. harmonic 
relative coordinates and velocities, x/j and v/j, and trans- 
form these expressions in terms of EOB coordinates and 
momenta. This transformation is facilitated by the fact 
that $e and $j being scalars, are actually expressed in 
terms of a basis of scalar combinations of x/j and v^. 
[Here x^ = — xJf, Vh, = dxh/dt = v'^ — v% are the 
relative, harmonic positions and velocities, considered in 
the center of mass system.] 
Let us use the notation 



x i = v h 



GM 



(3.16) 



and introduce X\, (A — 1,2,3) to refer collectively to 
these three scalars. The corresponding, natural EOB 
1,2,3 with 



scalars are X A , A 



X? = p 2 , xi=p 2 r . 



GM 1 



= - = u, 



(3.17) 



where, as above, p = p e /f- and q = r e /(GM). Note that 
all the scalars X A , X\ have the dimensions of a squared 
velocity. In other words X A /c 2 , X A /c 2 are dimensionless. 
In terms of this notation we have simply 



<T>' 1 
•i> E 



®%{X A ) 



(3.18) 

Therefore, starting from the known results for & E (X A ), 
$>j(X A ), it is enough to derive the transformation (taken 



(3.14) at a fixed, common dynamical time t h = t e ) 

X h A =f{X%) 



(3.19) 



to get the fluxes expressed in EOB variables. When 
PN-expanded the transformation (|3.19p has a polynomial 
structure, namely, 

X\ = £aB 1 X Bi + € 2 ^AB 1 B 2 X Bl X B2 

+e i UB 1 B 2 B 3 X Bl X B2 X B3+ 0(e 6 ) (3.20) 

Here e = 1/c is the PN expansion parameter and the 
structure of the 2PN-accurate expansion follows from the 
fact that Xa/c 2 ~ V 2 /c 2 is dimensionless. 

Actually, we have derived the transformation (|3.19[) 
by combining the two transformations that have been 
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explicitly worked out in the literature: (i) the transfor- 
mation between EOB (q e ,p e ) and ADM (q a ,p a ) phase- 
space variables [H, S 0] ; and (h) the transformation be- 
tween the ADM phase-space variables (q g . p a ) and the 
harmonic positions and velocities (q/^Vh) |4ll. [521. [53j . 

We give in appendix [E] the explicit forms of the various 
transformations (q e ,p e ) <H> (q a ,p Q ) «-» (<lh,Vh) we used, 
together with the explicit form of the resulting transfor- 
mation f|3 . 19[) . (|3.20p between the corresponding scalars. 



Cp, q (XA) = Ca 1 ...a p Xa 1 ■ ■ ■ Xa p + € 2 Ca 1 ... 



where we have q — p + 1 contributions, each one (using 
Einstein's summation convention) is a sum over all the 
indices A\ . . . A n it involves. Also the short-hand nota- 
tion 



XaXbXc 



X A 



BC. 



(3.22) 



will be adopted hereafter, when convenient. Note that in 
the multisummation 



.X, 



.x t 



the coefficient of (Ai)™ 1 (A 2 )™ 2 (A 3 )™ 3 (with n x 
n 3 = n) is 



S(ni, n 2 , n 3 ) C n 



22 



33 



(3.23) 
V n 2 + 

(3.24) 



where the symmetry factor S(ni, n 2 , n 3 ) is given by 

(nj +n 2 + n 3 )\ 



S(ni,n 2 ,n 3 ) 



nil n 2 \n 3 \ 



(3.25) 



In addition, as our basic variables are the EOB ones, we 
shall often, for brevity, suppress the index e (standing for 
EOB) on them: X A = X%. 

Before considering the higher PN corrections to the 
energy and angular momentum fluxes it is useful to re- 
call their leading order ("Newtonian order") expressions. 
They are easily deduced from the well known quadrupo- 
lar approximation (see e.g., 54]), namely 



_G_ 

5c? 



2G 



j(2) 7 (3) 



5 C 5 z ^ " ja ' 



with (in the center of mass) 



Iij = mixf l x{ + m 2 x 2 t x J 2 



<l~3> 



[IX X 



(3.26) 



(3.27) 
I 



By inserting the latter transformation in the results 
of Refs. [13 for the 2PN-accurate $'} we get the 
explicit expressions of $ E , in EOB variables. In order 
to better comprehend the structure of these results it is 
convenient to introduce a special notation for a general 
polynomial in X\. 

Given a collection of (symmetric) multi-index coeffi- 



cients C Al A 2 ... Ap 
< p < q and Ai 



■X A X Aw 



CaxA 



2 ■ ■ -ApAp 



i Ca x a 2 



1, 2, 3, we denote 
e 2 ^C Al ...A q X Al ...X Aq 



where 



(3.21) 



r 



using standard notation for symmetric and tracefree ten- 
sors. This yields 



G$ B 


1 




f GM 






~ <? 


5 ' 


\ r 




1 


V| 


(GM 










M 


= ? 


5 ' 


V r 



A 2 11 -2 

4ir r 

3 

2 0^.2 



J j I 1v l - 3f 



Note that these fluxes are both proportional to v 2 and 
contain a factor 1/c 5 (2.5PN order). In terms of a char- 
acteristic Newtonian velocity V (with GM/r ~ V 2 ,j = 
1/V), we have 



E 



^—T> 17 ~"— ■ 3 ' 28 

c 5 M c 5 



It will be also convenient to work with the quantities 

^5 



= — G$ E ~ vV 

V 

- c 5 <f> 1 ~ 
v M 



in 



(3.29) 



which have a finite limit when c — > 00 and in which one 
power of v has been factored out (so that they will be 
conveniently related to J- j [i). 

With this notation our 2PN-accurate results in EOB 
variables have the form 



^(X e ) = 

$5(x e )=j 



(GM 



CiAXa) 
3 



(3.30) 



Gi :3 {Xa) = Ca 1 Xa 1 



e Ca 1 a 2 Xa 1 Xa 
I 



where the explicit values of the coefficients entering 

£ 4 Ca 1 A2A3Xa 1 Xa 2 Xa3 (3.31) 



and Bi t3 (X A ) are listed in Appendix I A II and I A 21 Let us, for illustration, explicitly display here the leading or- 
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der contributions to $e and $,/ ("Newtonian order" 
namely 



C. Algorithm for decomposing &ej 



Finally, we need to compute the correspondingly 
rescaled version of the combined flux $fi °j B — — <fi e $j, 
Eq. (j2~l"2|) (with the EOB angular velocity e = <j> EOB ), 
in terms of EOB variables, namely 



e(Ncwt) 



4> 



:(Ncwt) _ 



8 [GM 
—v 

5 V r e 

8 (GM 
—v 

5 V r e 

8 (CM x : 

r{ — 



4p 2 -yp 2 
AX{ - yX 2 e 



j [ 2p 2 - 3p 2 + 2 



GM\ 



8 [GM 



j {2X1 - 3X| + 2X1) . 



c 5 G 



$| J = $| - GM(j> e $} . 



(3.32) 



Combining Hamilton equations for the angular motion, 
4> e = dH( co b)/dpct>, whose explicit expression is obtained 
from Eqs. f3TI3|) . i.e. 



-V = GAf<f = = 

eft f 2 KH {cS) 



J 



(3.33) 



with our above explicit expressions for <!>£ and 4>j, Eqs. 
p.30p . yields the following expression for 3>_ej, 



/ GM 



GM 



Ci, 3 (Xl) 



c 2 A / GM 



M-l( e S) V r e 



(3.34) 



where h has been defined in Eq. (|3.9[) . In Eq. (|3.34[) . and hence 
the factor {GM/r e ) 2 j 2 can be expressed in terms of the 
X%% Eq. (|3~T7| . since j = p /(GM) and 



GM 



i = -4 = (n e x p e )" 



= P 2 e~ (n e ' Pe) 5 



GM 



, ? 2 = X x e - XI 



(3.36) 



Similarly, one can replace the remaining factor 
(3.35) c 2 A/{hH {cS) ) = 1 + 0{e 2 ) in terms of the X^'s, namely 



„2 ±e 



c 2 A 



GMj hU 



1 + 



(off) 



X« + („_l)Xf + e 



3is 2 - v - 1 



(^33 _ ^13) 



+ ( i , + l)X 23 + -( l , 2 + z/+l)X 1 e 1 



(3.37) 



For instance, the leading order contribution ( "Newtonian order" ) to &%j reads 



ge(Newt) = Je(Newt) 
E J E 



GM \ 2 ge(Newt) 



5" 



GM 



GM 



4p 2 - ^rf r ) - (p 2 -p 2 ) ( 2p 2 - 3p 2 + 2 



8 /GM 



3 r 
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GM\ 



= — ^ 
5 



— ^ 

5 



GM 



GM 



Xf ( 4Xf - y X 2 e ) - {X{ - X 2 e )(2X 1 e - 3X 2 e + 2X1) 



-2(Xf) 2 - 3(X 2 e ) 2 + 5X^X1 + 2X^X1 - -XfXf 



-2X^ — 3Xf 2 + 5X^2 + 2Xf 3 — gX 23 



(3.38) 



Note that, while one could naturally factor u 4 = power of u = GM/r which one can naturally factor 



(GM/r) 4 in front of $ 



e(Ncwt) 
E 



it is only the third 



8 



out of $ 



e(Ncwt) 
E.J 



e(Ncwt) / 4 



E 



This difference is linked to the fact where 

(Newt) , 3 



[u was linear in Xa, while /u is 



that $ 

quadratic in the A^'s. 

When keeping the higher order PN corrections (which 
involve more powers of X A /c 2 ~ v 2 /c 2 ), the adimension- 
alized combined flux has the structure 



*!;j(X e ) = 



^4^ > 



(3.39) 



±{Xa) = Qa 1 a 2 Xa 1 Xa 2 + £ 2 Qa 1 a 2 a 3 Xa 1 Xa 2 Xa 3 + c^Q AiA 2 a 3 a 4 ,Xa 1 Xa 2 Xa 3 Xa 4i 



(3.40) 



the coefficients of which are listed in Appendix IA 31 

As indicated above, the first step of our new approach 
consists in separating out of $ e EJ either a factor Z\ = p 2 



or a factor Z 2 = rcW( oob )/<9r = — rp 7 
in terms of pi — Pi/fJ* and GM/r e = 
and Z2 respectively by 



As we are working 
l/q, we replace Z\ 



z,=f r = xi 

~ _ API 

Z-2 = —r e — - = 
dt 



dH( coh) 
or e 



4 I 



(3.41) 



which both have the dimensions of a squared velocity. In 
order to separate out a factor Z\ — Af or Z2 = Af from 
4>£j, Eq. (I3.39[) . we need to replace our basic set of scalar 
variables (A 4 , Af , Af) by the new set of scalar variables 
(Xf = ivi,Xf ,Xf = Z2). This is done by first expressing 
Xf = r e 9?i( eob )/9r e in terms of (XfjXfjXf) (by differ- 
entiating the EOB Hamiltonian (|3.6p with respect to the 
variable r e ) and then solving for Xf as a function of Xf , 
X§ and X 
have 



3 aiiu vi 4 . For instance, at the Newtonian order we 



(n 



(eob) 



Mc2) (Newt) 



!-2 

2 P 

1-2 



GM 
r 

IP% GM 
2 r 2 r 



so that 



Z. 



(Newt) 



X 



(Newt) 



3 



2q 2 



1 



<r q 
-X\ + X2 + X3 . 



(3.42) 



(3.43) 



Therefore, at the leading order, X\ can be solved in terms 
of X2 , X3 and X4 according to 



A'.,. = X-> A-; ! -.Y i +0(i 



(3.44) 



The extension of this result to 2PN accuracy is ob- 
tained by first computing Z 2 (Xi, X2, X3) to higher order, 
namely 



Z 2 = Xf(Xi,Xf,Xf) 
= -Ci,3(A"l) , 

where the coefficients of C*i 3 in 



(eob) 



dr e 



(3.45) 



Ci,3(X^) = Ca 1 Xa 1 + e 2 CA 1 A 2 XA 1 XA 2 + € A Ca 1 a 2 a 3 Xa 1 Xa 2 Xa 3 



(3.46) 



are listed in Appendix |B] Then one solves (perturbatively) for X\ in terms of X%, A3 and A4, starting with the 
Newtonian solution p. 441) . This yields 



X\ = XI + XI - X% + e z 2A| 3 + 3A; 



v — 5 



-XI 



v + l 



-X" 



1 + v 



-XI 



+e 4 [{2 - 6^)A 2 e 33 - 3(1/ - 3)A 3 e 33 + ^{v 2 + 7u - 63)A 3 e 34 + X -(y 2 - v + l)A 2 e 24 
+ 1(5^ + 8)A| 44 + IvXZu - \{v 2 + !/ + 3)A 2 e 34 - i(^ + hv + 1)A| 44 



(3.47) 



r 



where we have used the short-hand notation (already in- Here and below we find often convenient to use an ex- 
troduced in Eq. (|3.22[) for the variables Ai,A 2 , A3) plicit form for the polynomial expansion in powers of Xj's 

XfjK... = XjX e jX B K ... (I, J, K = 2, 3, 4) . (3.48) 
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(rather than a tensorial form CjXj + e 2 CijXjXj + . . . 
where one must take into account the symmetry factors 
associated with each term in the multi summations). 



Finally, by substituting the PN expansion of 
X£(X|,Xf,X£), Eq. dS17|), into the combined flux 
(I3.38|) . we get the expression of $£j in terms of Xf = 
(Xf,Xf,Xf). For example, at the Newtonian order, it 
suffices to replace Eq. (ET33|) into Eq. (pOS]) . This yields 



2X3X4 + -X2X3 



IX A — X0X1 



(3.49) 



r 



As anticipated, each term in this expression contains 
either a factor Z\ = X2 or Z2 = X4. It can there- 
fore be decomposed in the form (|2.15l) that we men- 
tioned above. Actually, there are many ways in which 
such a decomposition can be performed because the term 
— X2X4 = —Z1Z2 can be considered either as a part of 
$1^1 or of $2^2- 



We shall define the minimal decomposition (|2 . 15[) of a 
polynomial in the X/'s (which vanishes when X2 = = 
X 4 ) as the one of the form 



X2<f2(X 2 ,X 3 ,X i )+X i ^(X 3 ,X i 



(3.50) 



in which the coefficient of X^ does not contain any de- 
pendence on X2. (In other words, all the terms oc XJ 
are shuffled into the X2$2 contribution.) 

This minimal choice somewhat simplifies the expres- 
sion of $4, i.e., the coefficient denoted as $2 in Eq. 
(I2.I5[) - (I2.I7[) . In turn this simplifies both the radial com- 
ponent of radiation reaction and the Schott energy, be- 
cause, according to our above result (|2 . 1 T[) . these contain 
respectively dQ^/dt and $4. [Note the mnemonic rule 
that the indices get multiplied by a factor of two when 
passing from the notation of Sec. II to the notation here, 
Zi -> X 2 , Z 2 -> X 4 , $1 -> $2 and $ 2 -> $4-] 

For instance, at the Newtonian level, the minimal de- 
composition of $>ej reads 



5 \ T e 



X2 ( -X 3 — X 4 



- X 4 (2X 3 - 2X4) 



(3.51) 



while its 2PN-accurate generalization reads 

& EJ (Xi) = X 2 $ 2 (X 2 ,X 3 ,X 4 ) +X 4 $ 4 (X 3 ,X 4 ) 

(CM \ 3 r u 
— [X 2 $ 2 (X2,X 3 ,X4) + X4$4(^3,X 4 

where we found it convenient to factorize the term (GM/r e ) 3 in the above expression so that 

$2(X 2 ,X3,X 4 ) = (X|) 3 $ 2 , $ 4 (X 3 ,X 4 ) = (X 3 e ) 3 $ 4 



"4 , 



(3.52) 



(3.53) 



with 



+e 2 (^236 2 
+£ 1,105* * 24 



548 



76 



105 



24 21 

/1756 2 „ p 85494; 

v Xfo, H 

V 63 333 2 8 35 

159 ir 

"315 



5252 608 zuu 

105 33 105 23 105 

+24 2 80 

44 35 33 21 



44 1Q5 34 



^x ls + 

21 23 21 34 , 



105 

2835 333 105 223 



^ vX 22 3 - 



m Vx' 
if" 

4066 

— or — VJi.r 



!378 120268, 

_I/A 233 H 945" 233 



45 



32 ^x e 



1973 



Z^Xf 24 



45916 2 ^e 4066 

-3^5"^ A 234 35-^234 - 2X^^-244 - 

398 2 499 2 ,1496 14 597 25442 

+ ^T* * 244 " "63"^ * 224 + "siT^ 344 " ^5~^ 334 " ^15~* * 344 



10 



and 



892 



+ ■ 



105 
2428 

63 

1672 
945" 



z/ 2 X 



334 



668 
"35" 



289 
"35" 



701 
"35" 



945 V A 333 



I/ 3 X 



344 



"315^ * 334 



176 
105" 



V XnA A — 



857 
315 



16 



v^X 1 ' + —v 6 X. 



2:2:: 



^33 " 7T* ' X 



23 1 



16 



704 
105 
58 



is 2 X- 



33 



278 
105' 



256 
105" 



^ 2 X e 



568 
"35" 



45 



vX 



A ^x t 

35 4 

-^, 3 X 
35 



444 
3272 



1135 . 2 X e 

"63"^ ^344 



14597 
315 



v 2 X, 



94 5 333 



334 
1377 



34 105 

28f W 
" 315 334 " 



v z X. 



34 



538 
"35" 



9832 

"siT^ * 334 



vX- 



33 



94 5 333 35 



z/ 2 X 



333 



1363 
"315" 



-1/ XL.) — 



6536 
315 



vX 



3 14 



344 



(3.54) 



(3.55) 



It should be noted that $2 and $4 have the dimension of V 2 . Moreover, in the circular orbit limit X 2 = = X4 (for 
a later use) the above expressions reduce to 



32 

$2(0,^,0) = -^ 
10 

*4(0,Xf,0) = y*Xf 



e 2 X? 



e 2 X? 



28 



1313 
56 
269 
"21^ ~ "56" 



e 4 X 



/2291 



2195 213737\ 



33 



44 



+ e 4 X 



V 504 
/3041 



33 



V 1512 



168 
1377 
U2 l 



1512 J 
409 \ 
378 ) 



(3.56) 



r 



D. Minimal expressions of J> and -E( sc hott) 

Having obtained a particular, minimal decomposition 
of the 2PN-accurate combined flux $£,/(x, p) in the form 
(|2.15[) . namely Eq. (I3.52[) . we can now apply our general 
results (|2.17p . i.e., derive the corresponding minimal ex- 
pressions of T r and -E( sc hott)- Modulo the /x-rescaling 
(E = E/fi, p = p/fJ,), the prefactor (GM/r e ) 3 and 



($1 -> $ 2 , $2 -> $4), the second Eq. (|2.17|l yields the 
following minimal Schott energy per unit reduced mass, 

p(min) _ p(min) , 
-^schott — "schott/ P 



umi>) _ J_ / CM 

C 5 V r e 



schott 



p r $ 4 (xf,x 4 e ). 



(3.57) 



Note that the Newtonian order approximation to the 
(rescaled) Schott energy reads 



E. 



(min,Nowt) 1 16^ _ / GM \ 
schott =^—Pr{— 



1 16i/ 



GMY 
Te J 



O 



(3.58) 



where we used Eq. (|3.44[) to write the second form. The 
corresponding minimal expression of the (/i-scaled) radi- 
ation reaction is obtained from the first Eq. (I2.17[) . To 
write it explicitly, we first need to derive the value of the 
ratio pr/r. This is obtained from Hamilton's equation 

fo = £^b)_^ (Xjp)|irj (3 59) 



with 

where h is given by Eq. p. 91) above. The expression 
(I3.60[) for C is exact. Here, we shall work with its 2PN- 
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accurate expansion which is found to be 
C = l + e 2 C 1 , 2 (X%), 



imal (/i-scaled) radiation-reaction is given by 



(3.61) 



and the coefficients of C\$. are listed in Appendix iBl In 
terms of C, $2 and $4, the radial component of the min- 
I 



x"(cob) 
J~ r 



l.^MLp r ^ 2{ xi,xi,xi) 



d I ( GM 
dt { \~ 



(3.62) 



I 

Let us also recall that the azimuthal component of the and 
minimal (/i-scaled) radiation-reaction is simply given by 



( eob ) _ <T>( cob ) — 



j 



5^J 



1 / GM\ . 



jB lt3 (X%). (3.63) 



(Newt) _ P(j> , GM 



X 



(3.64) 



(where p<p is constant along the conservative dynamics), 
the time derivative of X3 and X4 are both proportional 
to r, e.g. 



(3.65) 



For illustration, let us display the leading-order ( "Newto- 
nian order" ) terms in these expressions. To get in explicit 

form the leading order expression of # eob) (x,p) we need 
to perform the time derivative in Eq. (|3.62[) by using the 
unperturbed (conservative) equations of motion. Here, 
we get some simplifications from having chosen $4 as a Re-expressing the result in terms of p r — r(l + 0(c~ 2 )) 
function of X3 and X4 only. Indeed, as Xf = GM j r e we get 



dxj Ncwt) = r_( 2 Pj>_ GA£ 
dt r \ r 2 r 



(X p)^-) = Ll„^Ml & ( 



21p 2 - 2lp 2 r 



GM\ 



O 



which, at this order, could alternatively be written in terms of velocities 

^2 



^(x,v)<»-I = ^<^)!,( 21 „ 



a -r) 



2 21f 2 - ) + O 



The corresponding, Newtonian order, results for J\/,(x, p)( Newt ) read 

3 - 



^(x,p)( Ncwt ) 
J- (x,v)( Newt ) 



1 8 (GM\ 

1 8 (GM) 2 ■ / 2 o n GM\ „ ( I 
■-=■— v± — <\> ( 2v 2 - 3f 2 + 2 l+O 

c° 15 



E 



(min) 
(schott) ^ 



1 . f GM\ 



(x, p) = ^p r I -— 1 (C A X% + e 2 C AB X AB + ^CabcXIbc) 
J- r (x, p) = l^p-p r (R A X% + e 2 R AB X AB + e 4 R ABC X ABC ) 



3 ~ 



- 1 (GM 

^ (X ' P) = ?( — / GM 



-{S A X C A + € 2 S A bX ab + e S A bcX abc ) , 



(3.66) 



(3.67) 



(3.68) 



The explicit 2PN-accurate versions of our minimal -E( S chatt)i j> an d T$ are given in Appendix [Cl and ID] They are 
expressed there in terms of X A = (Xf,Xf,Xf) and have the forms 



(3.69) 



where the coefficients C Al ... An , R Al ... A „, S Al ... An are explicitly displayed in Eqs. (|C2I) - (|D8[) . 
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The Schott energy as a function of X2, X3 and X4 (especially useful to study their limiting values along circular 
orbits) is given by Eq. p.57[) . while the radial and azimuthal components of the radiation-radiation force follow from 
Eqs. (EOH} and IpTM]) . i.e., 



Fr{X%,Xl,Xl) — -jr 



J0(X|, Xf, X\ 



PriTjXf + e 2 TjjXl T + e 4 T IJK Xf JK ) 



j(ViXf + e 2 V u Xjj + e 4 V IJK Xf JK ) 



(3.70) 



where / = 2,3,4 and the coefficients Tj 1 „j n , Vi x ...i n are explicitly displayed in Eqs. (ID12I) - (|D14[) . Note that if one 
wishes to express entirely in terms of X2, X3 and X4, the (rescaled) angular momentum term j should also be 
expanded in terms of X2, X3 and X4; the result is the following 



V X 1 ~ X 2 _ V X 3 ~ X t 



X'i 



X!, 



Wi 



xi - xt 



w 2 



( x § - x % 



(3.71) 



where 



W 2 



v+l 



XI 



44 



-X 



v+l 



31 



X 



4 " M ' 2 
4v 2 - 22v - 24 

X3444 



-X 



33 



v 2 + 8v + 1 



X\ 



32 

-3^ 2 + 4^ + 121 
32 

P- + 25v- 105 



X 



23 



4i/+l 



X 



32 

3(^ 2 - 1) 



3344 



16 

96v+ 16 
32 



X. 



3334 



8 



^■2333 2 2233 ' 



-X~2344 

16 32 
i^ 2 - 24z/ - 2 e v 2 + 3u + 5 

-^2334 + 



16 

3(^ 2 + i; 



XU 



X 



2244 



16 



^2234 



-12z/ + 27 



X 



3333 



(3.72) 



In the circular orbit limit these quantities reduce to 



vr 1 (o,x 3 e ,o) = -x 3 e 3, 



W 2 (0,Xf,0) 



-12i/+27 
8 



X 



3333 ■ 



(3.73) 



IV. NON MINIMAL CHOICES AND 
ASSOCIATED GAUGE FREEDOM 

Iyer and Will [H |H[ and later Gopakumar et al. HO] 
have shown that, at each order in the PN expansion, 
there is a multi-parameter arbitrariness in the construc- 
tion of a radiation-reaction force by the balance method, 
and that this arbitrariness is linked to the freedom in 
the choice of coordinate gauge. Let us briefly discuss 
how this arbitrariness enters our approach. First, it can 
be checked that our simplifying constraint (|2.9p that the 
Schott contribution to the angular momentum vanishes, 
•^(schott) = 0, corresponds to part of the freedom found 
by Iyer and Will. 

Indeed, one easily checks that within their approach, 
all the (non necessarily vanishing)parameters entering 
J( sc hott) are linearly independent, i.e., are unconstrained 
by the set of linear equations they obtained. Within 
our approach, this is immediately clear as we have ob- 
tained a solution with «/( sc hott) = 0, so that by choos- 
ing some given, general (nonzero) expression for J( sc hott) 
(such that J( sc hott) vanishes along circular motions) we 



will be able to straightforwardly construct a correspond- 
ing (minimal) radiation reaction force. [Indeed, the con- 
dition that ./(schott) vanishes along circular motion will 
introduce extra source terms in the equation (|2.11[) for 
J> and -E( sc hott) i linked to extra terms linear in Z\ and 
Z 2 in the right hand side of (|2.11[) , coming from an extra 
4>{8F,p) contribution to linked to SJ 7 ^ — — J( sc hott)-] 
This freedom in the choice of J( sc hott) is parametrized by: 
(i) one parameter (Ag) at the leading ("Newtonian") or- 
der, (ii) three parameters (AjjAgjA^) at the 1PN order, 
and (iii) six parameters (A^, A33, A44, A23, A24, A34) at the 
2PN order. The general form of J( scho tt) = ^(schott) //■* can 
be written as 

?(non min) 1 ~ ~ / GM \ , j 
J ( S chott) = ^PrP4> [j- ) ■ ( A 

+e 2 A/X / + e^jjXjXj) , (4.1) 

where the free gauge parameters parametrize the coeffi- 
cients of a general polynomial in Xj = (X2, X3, X4). 

Note that these parameters were indicated differently 
m previous papers d |H, Hi]. In particular, the single 
J-related parameter Aq at leading ("Newtonian") order 
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was previously notated as 



GII 



JW 



-BD 



(4.2) 



and was normalized so that Aq = (8/5)^/3^ /7 . 

Besides the parameters associated with the (non min- 
imal) choice of a non vanishing t/( sc hott)i there are fur- 
ther arbitrary parameters which, in our approach, cor- 
respond to further non minimal choices in the construc- 
tion of -E(schott)- Indeed, our general result (|2.17p shows 
that the arbitrariness in the coefficient $2 of Z 2 in the 
decomposition (|2.15|) . will directly affect £*( schott), an d 
then T r . [Given a choice of $2, compatible with (|2.15l) . 
the corresponding $1 is uniquely determined.] As dis- 
cussed above, the arbitrariness in $2 is parametrized by 
a general term oc Z\ = p 2 = A" 2 . In terms of the rele- 
vant basis X 2 , X3 , X4 (with Xi oc Z\ and X4 oc Zi ) the 
arbitrariness in $4 ~ $2 in Eq. (I3.41[) is of the general 
form 



(non min) 



{Xj) = vX 2 (A 7 + e 2 XjXj 
+e i \ J IJ X I X,j) , 



(4.3) 



to -E( S chott) of the form 



corresponding to an additional non minimal contribution 

(4.4) 



E 



(non min 
(schott) 



-3 



C \ TV 

r 2\E v 1 A\E 



-e z \fX I + e i \fjX I Xj) 



This expression shows that the additional gauge-freedom 
associated with such non minimal choices in the Schott 
energy is parametrized by: (i) one parameter Xq at 
the leading (Newtonian) order, (ii) three parameters 
(Xf , Xf , Xf) at the 1PN order, and (iii) six parameters 
(Af 2 , Af 3 , Af 4 , Xg, Af 4 , Af 4 ) at the 2PN order. 

In terms of the notation of [40| (if we approximately 
identify their Lagrangian framework with our Hamilto- 
nian one) these parameters correspond, respectively, to: 
(i) "3, (h) £2161, & and (iii) V2 , V>4 , V>6 , V>7 , V's , V>9 , i-e. to 
the following contributions (oc r 3 ) to the Schott energy 
considered in 14011: 



E 



(non min, Newt) 
(schott) 



16 G 2 M 2 



E 



(non min,lPN) 
(schott) 



16 G 2 M 2 



r 3 (&v 2 +^r 2 



Hi 



GM 



^(nonmin^PN) 16 G 2 M . 3 , 4 , , 2-2 
C T 

j GM , . 4 , . 2 GM 
+ipev 1- W7 r + V8f 



(schott) 



. (GM\ 



(4.5) 



Summarizing: the arbitrariness in the construction of a 
radiation-reaction force is parametrized by the parame- 
ters Aq, Aj , Xj j . . . entering the (non-minimal) Schott 
angular momentum (|4.1[) . together with the parame- 
ters Aq 3 , A 7 ?, Xf il2 . . . entering the (non-minimal) 0(p 3 ) 
Schott energy (]4.4[) (expressed as a function of X2 = p 2 , 
X3 and X4). It is easy to see that the number of arbitrary 
parameters entering the nPN order is equal to 



n + 2 
n 



n + 2\ (//+ Di// +2) 



(4.6) 



for each one of these sources, with ao = 1, ai = 3, a-i = 6, 
03 = 10, etc. 



V. SOME APPLICATIONS OF OUR RESULTS 

A. Schott energy along quasi-circular inspirals 

Recently, Damour, Nagar, Pollney and Reisswig [55[ 
have compared several different functional relations E{ J) 
between the energy E and the angular momentum J of a 
binary system evolving along a radiation-reaction driven 
sequence of quasi-circular orbits. In particular, they com- 
pared a relation i? NR (J) obtained from accurate numer- 
ical relativity (NR) simulations, to several of the rela- 
tions E EOB (J) that can be derived from EOB theory (un- 
der various approximations). Actually, the NR relation 
,E NR (J) computed in Ref. [55| was obtained by defining 
the NR energy i£ NR and the NR angular momentum as 
being their initial values minus the time integral of their 
respective NR fluxes, and $j R (as recorded at in- 
finity). In view of our general balance equations (|2.7p . 
we see that (modulo numerical errors) the NR energies 
and angular momenta can be identified with the sum of 
the system plus Schott contributions: 



E NR (t) = £ (system) (x(i) )P (i)) +£ (Bdrott) (x(t),p(t)) 

J NR (i) = ./(system) (X(*),P(*)) + J( S chott)(x(i),p(i)) . (5.1) 



On the other hand, one of the tenets of the current im- any moment, to minus the angular momentum flux $j. 
plementation of the EOB formalism is to require that the 
^-component of the radiation-reaction force be equal, at 
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In view of the second Eq. (|2.8p . this means that the 
EOB formalism has chosen a "gauge" where 



«tt)(xW, P W) = o. 



(5.2) 



In view of this, it is consistent to identify the instanta- 
neous NR angular momentum J NR (t) with the EOB one 
J EOB , which indeed measures the angular momentum of 
the system, J (syst 

cm) } • 



J«K (t) = J*U« (x(t)i p( t) ). 



(5.3) 



By contrast, in view of the first equation (15.11) . the EOB 
measure of the total energy of the system, defined as 

B EOB (x(t),p(i)) = H (eob) (x(t),p(t)) - Mc 2 

= E (syatem) (x(t),p(t)) , (5.4) 

cannot be simply identified with the NR computed en- 
ergy _E NR . Indeed, one expects the relation 

E NR (t) = ^OB (x(t))pW) + £ EOB tt)(x(t))p(t)) (55) 

In conclusion, as was already pointed out in Ref. (Hlj . 
the NR-derived functional relation i? NR (J) should dif- 
fer from the EOB derived one E EOB (J) by the quan- 
tity -E( sc hott) (*) > re-expressed in terms of the correspond- 
ing instantaneous angular momentum J(t) = J NR (i) = 
J EOB (t). 

Our results provide, for the first time, the explicit an- 
alytical value of -E( sc hott)j namely the first of Eqs. ()3.69p 

I 



(see Appendix [C]). Note that -E( sc hott) is proportional 
to p r , which stays rather small all along the radiation- 
reaction driven sequence of quasi-circular inspiralling or- 
bits, including most of the subsequent plunge phase (see 
Fig. 1 of [2J). The smallness of p r further implies that 
the numerical value of S( sc hott) is approximately gauge- 
invariant during the inspiral and the plunge. Indeed, Eq. 
above shows that the general non-minimal contri- 



bution to i?( sc hott) contains an overall factor p^, instead 



of the corresponding factor p r in E 



(min) 



The ratio 



(schott) ' 

^sIoTt) n) / E (Tchott) is therefore generally expected to be 
numerically of order p 2 , and hence small during the in- 
spiral (and even the plunge). 

In addition, during the inspiral, i.e., before crossing 
the Last Stable (circular) Orbit (LSO), the dimension- 
less scalar X%/c 2 = Z 2 /c 2 = (r e /c 2 )dl-L( eo b)/dr e will also 
be numerically small. [Indeed, the orbital radius r e (t) 
approximately stays at the bottom of the effective radial 
potential "H( e ob)(j) r e) during the inspiral.] The numeri- 
cal value of -E(schott) during the inspiral can then be ap- 
proximately evaluated by neglecting X\ in $4 (X| , X|) . 
This leads to an approximate, simplified expression for 



-E'(schott) . 

radius r, 



along the inspiral, as a function of the EOB 



^(inspiral) 
(schott) 

i.e., explicitly 



^Kt)- ? jir*4(*!.0), (5, 



fp(inspiral) 
(s, 

1 



3 r 



3024 



(6082^ + 37179^ + 3272) 



1-7^(807 + 352^) , 
168 \ r, 



GM 



(5.7) 



Note that -Evgchott) (*) * s ne 9 a ^ ve (because p r ~ r < 
during the inspiral). It would be interesting to take into 
account the modifications of the EOB /NR comparison of 
Ref. [55[ introduced by the presence of the Schott contri- 
bution to the energy (especially during the late inspiral 
and the plunge). This might allow one to refine the con- 
clusions of Ref. [55j and to extract some information 
about the exact form of the EOB Hamiltonian. 



fied form 

T r = (5.8) 
= / ■ (5.9) 

For instance, at the Newtonian order they argued that 
the choice 



«BD = QTW = /?2GII 



10 



(5.10) 



B. About the radial component of 
radiation-reaction 

When Buonanno and Damour [2] incorporated 
radiation-reaction effects in the EOB formalism, they 
suggested that it is possible to use the radiative gauge 
freedom to put the radiation-reaction force in the simpli- 



of one of the two free gauge parameters entering jr( Newt ) 
ensured the vanishing of the radial component J r r Nowt ' ) . 
This statement is correct. However, this specific choice 
of «bd = otxw — P2G11 conflicts with the second require- 
ment (|5.9[) that J-0 be identified with minus the angular 
momentum flux. Indeed, our results above (as well as 
the previous results of Iyer and Will) show that the sim- 
plifying requirement (|5.9p actually determines the value 
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of half of the free gauge parameters entering Ti- More 
precisely, they determine the values of the parameters 

( n = 0. !' 2 ) enterin S J ( s n c hou) in) , Eq. 6U (namely 
Xj j = 0). One the other hand, as pointed out in 
Sec. IV above, the Newtonian order J( sc hott) -related pa- 
rameter Xq happens to be proportional to the parameter 
<5bd = a\w = @2Gii which needed to take the nonzero 
value (|5.10p . We see therefore that the choice (|5.10|) cor- 
responds to a non-minimal (i.e., non vanishing) value for 
J( sc hott)) in conflict with the second, simplifying require- 
ment (|5T5|t . 



In view of this result, we henceforth advocate to in- 
corporate radiation-reaction in the EOB formalism by 

I 



consistently enforcing the minimal choice 

J> = -*j, (5.H) 

corresponding to 

J ( i° h B ott) (x(t),p(t))=0, (5.12) 

i.e., Xj j = 0. This choice necessarily implies a nonzero 
value for T r . In particular, if we also require the second 
minimal choice, 

if s ° B ott) (x(i),p(t)) = if ° h B ot ?f (x(i),p(i)) , (5-13) 

we have seen above that T r is completely determined, 
and has the form 



?r — -^Pr- — ~—(RaX% + e2 RAsX\ B + ( 4 RabcX abc ) 
c° ri 



(5.14) 



I 

where the coefficients Rai...a k are listed in Appendix ID] comes 
If we consider the case of a quasi-circular inspiral we can 
neglect X| = Z\ — pf,, and replace Xf = p 2 by the 
expression obtained by setting Xf and XJ to zero in the 
relation (|3~47l) . 

Specialized along circular orbits, relation (|3.47j) be- 

I 



X[ cm) = XI + 3e 2 X£ 3 - 3(v - 3)e 4 X 3 e 33 . (5.15) 
This leads to the following approximated form of T r 



f-Onspixal) „ "32 (GM) 3 

~ c 5 3 Pr r* 



1 - 2 lo (1133 + 94 ^ 



GM 



1 



15120 



(-175549 + 322623^ + 70794^ 2 ) 



GMY 



(5.16) 



It might be useful to record the value of the ratio between the inspiral value of T$ [obtained by replacing Xf — > 
T r and T§ during inspiral. To this end, we first note that and XJ — > in Eq. (|3.63|) ] reads 



r. 



r n::M.i,;Hi _ I' 32 ( ( r M ^ 



so that we have the ratio 

^.(inspiral) _ - 

= -gmOl 



T 



(inspiral) 



c 5 5 V r, 



7/2 



1 -— (588^+ 1247) 



^^(-89422 + 153369^ + 9072^ 2 ) \ > 1 



227 1957\ ( GM\ 2 /753 2 165703 25672541\ / 'GM 



(5.17) 



1+ V 140^ + 1680 J V r e J 6 + \m0 U + 70560 V 5080320 J \ r, 



(5.18) 



This result is consistent with Eqs. (3.14), (3.18) of 2] with the value o,bd = (i.e., Aq = 0). We leave 
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to future work a detailed study of the consequences of 
incorporating in the EOB formalism the non-vanishing 
value of T r advocated here. The preliminary comparison 
performed at the end of Sec. V in Ref. @ (between 
using Jv/J 7 ^ = and J>/J-0 = f/(r 2 <fi)) indicates that 
the effect of the more consistent value of J>/J-0 found 



here will be small. However, modern use of EOB theory 
aims at a very high accuracy in the phasing, for which 
the new value of J- r will probably have a significant effect. 
Let us also recall that along circular orbits, one finds (at 
2PN order), using J^ circ) = p|/r 2 and Eq. (j5TT5j) . 



p = ^GMr e 

and hence, for fi( clrc ) = dU( coh) / dp,p\ c - m , 

GMO( circ ) _ (GM 
c 3 \c 2 r e 



, 3 GM 2 3,, 



(5.19) 



.3/2 



vGM 2 3 
2T7 e + 8 ( ^ 5) 



(5.20) 



The latter equation implies the following expression for the dimensionless frequency parameter x, i.e. 



/GMfi( circ ) 

V C 3 



2/3 



- f GM \ 



V c2r e 7 



i/GM 
3 c 2 r P 



36^ 



"(SO' 



inverting (perturbatively) this relation yields 



GM 



= x vx 

3 



—vx 
4 



(5.21) 



(5.22) 



so that, in terms of x we have 



^S(inspiral) 
(schott) 

x-(inspiral) 
J~ r 

-£-(inspiral) 



1G 



-vx 3 p r 



65 



32 c 3 

Tgm 



4 ~ 

vx p r 



1 + 



269 
21^" ~56~ 
494 



/7769 

1 

V 1512 



7543 
336 1 



32 2 7/2 
— C VX ' 

5 



1 



105 
35 



1133\ 
280 ) 
1247 \ 

336~ r 



/3071 2 



409 \ 
+ 378/ 
55577 



V 280 
65 
18* 



1680 



175549\ 
15120 J 



9271 
"50?' 



44711 
9072 



(5.23) 



The latter expression of jK lns P iral ) as a f unc tion of the frequency parameter x agrees with well-known previously 
derived results (see, e.g., Eq. (4.18) in [Hj]). 

C. Hyperbolic orbits: conservative aspects 



Up to now, the EOB formalism has been applied 
only to the description of radiation- reaction driven quasi- 
circular orbits, because these are the orbits of greatest 
relevance for the current network of ground based grav- 
itational wave detectors. However, we anticipate that it 
will be useful to apply the EOB approach to other or- 
bits, such as elliptic orbits, but also hyperbolic ones. It 
is now possible to do so because we have provided above 
a description of radiation-reaction along general motions. 
Here, we shall consider the case of hyperbolic motions, 
and focus on the effect of radiation-reaction on the an- 
gle of scattering of a gravitationally interacting binary 
system (considered in the center of mass system) . 

Before taking into account the additional effects of the 



I 

radiation-reaction force Ti, let us consider the conser- 
vative dynamics of hyperbolic encounters (at the 2PN 
accuracy). We recall that, at the 2PN accuracy, the rela- 
tive motion in the orbital plane, r(t), (f>(t) is described (in 
any PN g auge; h armonic, ADM or EOB) by equations of 
the form [M&tEl 



dr 



-s =2E' + - 
dt J i 

2-^3 



2' u 2 y 



2 



r 
■ Ri 



ft** = j' 

di 



,Gi 



L Go 



4^; 



1 + e — + t— + e 1 — 



(5.24) 
(5.25) 



Here we have used the scaled variables (f = rj (GM)), t = 
t/(GM)), and the prime on any quantity Q denotes a 
multiplicative modification by higher PN terms of the 
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type Q' = Q(l + qi€ 2 + q2e 4 ), where qie 2 , q2t A (as well as 
the coefficients R p e q , G p e q above) are polynomials (with 
^-dependent coefficients) in the dimensionless quantities 
E/c 2 and l/(jc) 2 . For instance, at the 1PN accuracy, 
and in harmonic coordinates 59] 



2E' = 2E l + |(3^-l)^ + ofl N 



2' = 2 i + ( 7 „-6)* +0 (£ 



-0 



f = j l + (3„-l)§ + (£ 



(5.26) 



Note that (j 2 )' is not the square of f. Many previous 
investigations [IH [57l - l59| were interested in describing 
the motion as a function of time. Here, we shall instead 
focus on the shape of the orbit, i.e., f as a function of <f>. 
This is obtained by eliminating dt between Eqs. (|5.24[) 
and (|5.25p . Introducing the dimensionless variable 



(J'f 



(5.27) 



leads to a first-order differential equation for u((f>) of the 
form 



^ ) = 2E'{j') 2 + 2'u - l'u 2 + e 2 U 3 u 3 + e 4 U 4 u 4 + e 4 U 5 u 5 



(5.28) 



where all coefficients (£', {j') 2 , 2', 1', e 2 J7 3 , e 4 f7 4 , 
e 4 £/s) are dimensionless. One can then reduce the 
above equation to a Newtonian-looking equation by a 
suitable change of (inverse) radial coordinate. Indeed, 
by appropriately choosing the dimensionless coefficients 



e z Ci,e^c 2 ,e^c 3 m 



(5.29) 



one can get (modulo 3PN terms) an equation for u((f>) of 
the simple form 



^ ) = 2{Ej 2 )" + 2"u - l"u 2 , 



(5.30) 



where the double prime indicates further multiplicative 
modifications by higher-PN terms of the usual coefficients 
entering the Newtonian-order equation for vS N ^ = j 2 /r, 
namely 

(^^j 2 = 2Ej 2 + 2uW - (/>) 2 . (5.31) 

The general solution of the latter (Newtonian-order) 
equation is the well known polar equation of a conic, 



(</») = l + e (Ar) cos0 



(5.32) 



with eS N ~) = J 1 + 2Ej 2 . By contrast, the general solu- 
tion of the modified Eq. (|5.30[) will be of the form 



(5.33) 



u{cj)) = C I 1 + e cos 



K 



I 

where (j>o is an arbitrary integration constant and where 
e, C and K are functions oiEj 2 , E/c 2 and l/(cj) 2 which, 

respectively, reduce to \Jl + 2Ej 2 , 1 and 1 when 1 /c 2 — > 
0. Note that the quantity K measures the periastron 
advance 



K= — 
2tt 



1 



(5.34) 



where $ denotes the period of <f> (i.e., u(<j) + $) = u(<p) 
in the elliptic case; see below the definition of $ in the 
hyperbolic case), and where k is the usual notation for 
the relativistic contribution to periastron advance. It is 
given at 2PN by [57| 



HE, 3) 



(cj) 2 
35 

T 



ijcf 



E 
c z 



O 



(5.35) 



[See Ref. [60] for the 3PN accurate value of k] . Here, we 
work with the analytic continuation (in E) of the function 
k(E,j) from the elliptic-like case (where E < 0) to the 
hyperbolic-like one (E > 0). Note that we can further 
simplify the result (|5.33|) by modifying the leading-order 
coefficient 1 in the parenthesis appearing on the right 
hand side of Eq. (|5.29|) so as to absorb the coefficient 
C = 1 + 0(c~ 2 ) in a rescaling of u. In other words, there 
exist coefficients 1' = l + 0(£Vc 2 ) + 0(l/(jc) 2 ), and e 2 ci, 
e 4 C2, e 4 C3 such that the polar equation f((j>) of the orbit 
takes (at 2PN order) the form 
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£ i' 



i 2 

2- J 

e ci — 

r 



e 4 c 2 



T =1 



e cos ■ 



A 



(5.36) 



This form is valid in any PN gauge (harmonic, ADM 
or EOB). We will give below the explicit values of its 
coefficients in the EOB case. In this form the two coeffi- 
cients, e and K entering the rhs acquire a gauge-invariant 
meaning. This is well known for the periastron advance 
parameter K (when it is considered for the elliptic-like 
case), but this is also true (when considering asymptot- 
ically flat gauges) for the "eccentricity" e (when consid- 
ering the hyperbolic- like case). Indeed, when considering 
hyperbolic orbits the lhs will vanish both in the infinite 
past (incoming state, f — ¥ +oo) and in the infinite future 
(outgoing state, f — > +oo) so that (choosing the integra- 
tion constant (j>o = 0; location of the periastron) <fi will 
evolve from in the infinite past to <f> + in the infinite 
future, where 4>-(= — 0+) and 4>+ are the two solutions 
of 



Both the scattering angle x an( i the periastron pre- 
cession parameter K are gauge-invariant physical quan- 
tities that can be expressed as functions of the two basic 
gauge- invariant dynamical parameters E and j . We see 
therefore from Eq. (|5.40l) that e can also be considered 
as a gauge-invariant quantity, and can be, in principle, 
expressed as a function of E and j 2 . [We shall give below 
some explicit integral definitions of the functions x(i?, j) 
and K{E,j) from EOB theory] 



At the 1PN accuracy, the invariant eccentricity e coin- 
cides with the eccentricity denoted as eg in [59[ (see Eq. 
(5.7) there, which is of the form (|5.36jl ). The expression 
of e 2 in terms of E and j 2 is given by (see Eq. (4.13) in 



1 + e cos — = , 
A 



(5.37) 



i.e. (we are in the hyperbolic case where e > 1) 

<t>± = iAarccos (~^) ■ ( 5 - 38 ) 

The (center of mass) scattering angle, x (taken with a 
positive sign) is related to <p± via 

X + 7T = 0+-0_ = A0 (5.39) 

so that we can write x m terms of K and e according to 

r 



X + 7T = A(f) = 2Aarccos 



(5.40) 



J 



-2 2 i 

e = e„ = 1 



2Er 



1 - 



6 



15\ E 



(5.41) 



We have determined the extension of this relation to the 
2PN accuracy by using results in the literature on the 
"quasi-Keplerian" parametrization of the 2PN motion 
[571 158| , namely 



nit — to) — u — e f sinu H r sinu + —riv — u) 

c 4 c 4 



= a r (l — e r cos u) 



K 



—r sin 2v H — r sin 3v 

c 4 c 4 



(5.42) 



where 



v = 2arctan 



l + e 4 
1 — e-j 



1/2 



U 

tan — 
2 



(5.43) 



Here the "eccentric anomaly" u (and its analytic contin- 
uation u mentioned below) should not be confused with 
the gravitational potential variables u — GM/r e , u, used 
above. 

The form written here corresponds to elliptic-like mo- 



tions (E < 0). However, similarly to the Newtonian 
case (which is recalled in Appendix [FT) the corresponding 
parametrization of hyperbolic-like motion is obtained by 
the simple analytical continuation 



iu . 



(5.44) 



which accompanies the continuation of E from negative 
to positive values, as well as the continuation of the var- 
ious eccentricities et, e r , from < 1 to e, > 1. In 
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addition, n 2 ~ GMja z r and a r ~ —GM/(2E) are contin- 
ued from positive to negative values. In this continuation 
the angular variable v remains real. The radial motion 
equation becomes r = a r (l — e r coshu), so that the out- 
going and incoming states are described by u — )■ ±00. 

This corresponds to finite values of the real angle v 
given by 



v± = ±2arctan 



so that (choosing O = as above) 



1 



1/2 



(5.45) 



where 



B 2 



1 - 6 



-|(10z,-17)A. 



(5.48) 



This leads to several possible ways of computing the 
scattering angle \ as a function of E and j 2 , at the 2PN- 
accuracy. A first form would be obtained from Eq. (|5.40p 
without any re-expansion, i.e., 



<t>± . U ■ o , 50 • o 

— — = v± H — j sin 2v± + —r sin 6v± . 
K c 1 c 1 



(5.46) 



Taking the cosine of this result, and using the 2PN ac- 
curate expressions of e^, f$, as functions of E and 
j 2 [58| then leads to the following explicit 2PN-accurate 
expression for e(E,j) 

e 2 (E, j) = l + 2Ej 2 [1 + e 2 B 1 + e i B 2 + 0(e 6 )} (5.47) 



x(E 1 j)+n = 2K(E,j)&rcos - 



1 



(5.49) 



where K(E,j) = 1 + k(E,j) is written in Eq. ([535)) 
above, and e 2 {E, j) is the polynomial in E and j 2 written 
above. 

Alternatively, one might consider re-expanding the re- 
sult (|5.80p as a straightforward expansion in 1/c 2 . This 
leads to 



7;X(E,j) = arctan 



'2Ef 



e 2 A x + e 4 A 2 , 



(5.50) 



where 



3 - \/2£?' 2 

Ai = ^<t>° + (E,j) i = \(v - lb)Ej 2 - 12] 

j J ' Aj 2 (l + 2Ej 2 ) lK ' J 

A 2 = <jP + {E,j)A 2a + A 2b , 



(5.51) 



and 



A 2a — —k 
3 



A 2b = 



(5-2!/) - 5(7-2!/) 



2£j 2 



2(3i/* + 30i/ + 35)E 3 f + [y l - 838^ + 2593).E z j 



2-4 



32(l + 2£j 2 ) 2 j 4 
-32(28i/ - 95)iy 2 - 240i/ + 840 



r 



(5.52) 



6° (jE, j) = arccos — 



1 



1 + 2Ef 



(5.53) 



One can also consider the PN expansion of tan ¥ , namely 



tan 



X(E,J) 



2Ef 



1 + 2Ef 



2Ef 



e 2 A 1 



A 2 + 



,A\ 



2Ef 



(5.54) 
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Beware that the straightforward PN expansions of 
k(E,j) and x(E,j) are badly convergent because of 
the presence of a singularity (where k(E,j) — > oo and 
x(E,j) — > oo) along the sequence of unstable circular or- 
bits. Let us recall that, in the (E, j) plane the sequence 
of circular orbits is defined by parametric equations of 
the type (when v — > 0) 



E(x) 



Cj(x) 



l-2x 
7! -3a; 



l + 0{y) 



\/ x{\ — Zx) 



0(y). 



The orbits we consider here (either elliptic-like or 
hyperbolic-like) lie between the two branches defined by 
the parametric equations above: the lower branch of sta- 
ble circular orbits (corresponding to < x < xlso( z/ )j 
with xlsoM = ^ + 0(v)), and the upper branch of 
unstable circular orbits (xlsoM < x < :elr(^), with 
2 ; lr(z / ) = 3 + 0(v)). [Both branches meet at a cusp cor- 
responding to the LSO.] Many of the functions of E and 
j that we consider here (and notably k(E,j) and x(E, j)) 
become singular on the upper branch. It might then be 
better to work with the PN expansions of related vari- 
ables that are regular on the upper branch, e.g. related 
functions that smoothly vanish there instead of blowing 
up. [When considering the zero-eccentricity limit, this 
strategy was used in Refs. [3J, |60j, which replaced the 
singular function K(E,j) by the smoothly vanishing (57j 
function K- A (E,j)] 

Let us finally note that the EOB formalism gives an 
exact integral form for the scattering angle. Indeed, ap- 
plying the Hamilton- Jacobi method to the EOB Hamil- 
tonian leads to a separated action of the type 

S(eoh)(t, r , <t>; E,p^) = -Et + 



drp r (r; E,p^) , (5.55) 



where p r (r; E^p^) is obtained by solving the equation 
T~L{cob) = E, or, in terms of the /x-reduced effective energy 
^(cff) (using also f = r/(GM)) 



c 2 f 2 c 2 B(f) 



This yields 



(5.56) 



(5.57) 



The orbit 4>{f) is then obtained from using dS^o^/dp^ = 
4>a =constant. Setting 0o = yields 

0(f) = -(GM)J- j dfp r =±J dfn(f;j,H {cS) ) , 

(5.58) 



where 



" '■; 
T3 - 



-A{r) (l + &) 



(5.59) 



It is useful to re- write this result in terms of the inverse 
radius u — 1/f — GM/r. 
Introducing 



^A(u)B{u) 



H 



%a _ A(u) (c 2 + fu 2 ) 



we have 



u)=± duJ7(u;j,H (eff) ) 



(5.60) 



(5.61) 



The function U(u) is defined as a real function in the 
classical domain where the function appearing under the 
square root in its denominator, say T)(u\ j, T-L^q)) = T>(u) 

V(u) = ^ - A{u) (c 2 + jV) (5.62) 

is positive. In the elliptic-like case {E < 0) this is the 
case in an interval of the form < u min (E,j) < u < 
j), where w m i n and u max are two positive roots 
of V(u). In the Newtonian approximation 2?(u)( Ncwt ) = 
2E + 2u — j 2 u 2 , these two positive roots are 



(Newt) 



(E,j) 



l-Jl + 2Ef 



l + x /l + 2Er 

J 2 



(5.63) 



Then the angular period $ = 2ttK is given by an integral 
over the interval [u m i n ,u max ], namely 



irK(E,j) 



<HEJ) 
2 



duU(u;j,n {eS) ). (5.64) 



When one continues E from negative to positive values, 
the analytic continuation of u m i n (E 7 j) stays real, but be- 
comes negative. However, nothing wrong happens to the 
integrand, and one can still consider that the real integral 
above defines K(E,j) in the hyperbolic-like case (E > 0, 
i.e., "H( c ff)/c 2 > 1). [In terms of the usual radial variable 
f = 1/u this means that one is taking an integral which 
goes beyond r = +oo to formally extend to negative val- 
ues of the variable f] 

By contrast, the scattering angle \ is directly defined 
in the hyperbolic-like case (E > 0) by an integral over 
the interval < u < u max (£', j), namely 

x{E,j) 7T _ A<t>(E,j) 
2 2 2 
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1 



(5.65) 



Here the interval < u < u max (E,j) corresponds to the 
radial interval f min < f < +oo, where f min = l/w max 
is the minimum of f (periastron). By comparing Eq. 
(I5.65|) with Eq. (|5.64l) we see that while K is given by 
a complete integral (i.e., a period integral, between two 
successive roots of U(u)), \ is given by an incomplete 
version of the complete integral (going between a root 
and u = 0, which is an intermediate point). This explains 
why the PN expansion of x(E,j) nas a more complicated 
analytical structure as a function of E and j 2 [involving 



arctan(l/y2£j 2 )], than K(E,j). 

Let us finally indicate how one can rather easily com- 
pute the explicit quasi-conical equation (see Eq. (I5.36p ) 
of the orbit in EOB coordinates. Let us consider the 
squared differential of the polar angle, dtp 2 = U 2 (u) du 2 . 
We wish to transform it, by a (2PN-accurate) change of 
u variable of the form 

u = u + e 2 au 2 + e 4 bU 3 + ^cU 4 + 0(e 6 ) , (5.66) 

so that it simplifies (modulo 0(e 6 )) to a form involving 
a quadratic polynomial in u as denominator, i.e. 



f/ =J 2 



D(u) (du) 2 



J' 



(c 2 + e) - A(u)(c 2 + j 2 u 2 ) 
(du) 2 

'2RTT 2 ' 



£ + 2(111- j 2 j3 u 



(5.67) 



Here, D(u) = A(u) B(u) — 1 — e 2 6vu, and we introduced 
the new energy measure e (not to be confused with the 
PN ordering parameter e = 1/c) 




= 2 E+— E z ) + — [E + — E 



It is easy to check that the choice of coefficients 



■ = -l-4}b(17-10«0 + o(i 

b = |(i + a/) + o(l 



(5.68) 



c = + O 



1 



(5.69) 



in Eq. (|5.66[) does yield the simple u-form indicated in 
the second Eq. (|5.67|) . The coefficients a and /? entering 
the quadratic u-denominator e + 2 a u — j 2 f3u 2 are then 
found to be (at 2PN accuracy) the following functions of 
E and j: 



a(E,j) = l-^e + o(± 



1 



K 2 (l + k) 2 



(5.70) 



[The latter result for (3, that we explicitly checked at 
2PN, must hold to all PN orders.] 

The integration of Eq. (|5.67p then yields 



u = (u) ( 1 + e cos — 
K 



where, denoting by u\ and u-i (u\ < U2) the two roots of 
the quadratic u-denominator, 

e + 2au - f/3u 2 = j 2 (i(u - ui)(u 2 - u) , (5.71) 

we have 

Ui + u 2 _ U2 — U\ 
(") = 7; ' e 



Ul + U 2 



This yields 



and 



a aK 2 



e 2 -l=ef 



3 2 P J 2 
f> ej 2 



m 2 K 2 a 2 



(5.72) 



When inserting in Eq. (|5.72p . the expressions of e 
(Eq. ([5T681) ). K (Eq (|535)» . and a (Eq. ([570]) with the 
first equation (|5.69[l ) . one finds, after PN reexpanding 
e 2 (E,j) the same result as Eq. (|5.47|) above. 



D. Hyperbolic orbits: radiative effects 

Having explained the various ways in which one can 
compute the scattering angle x as a function of E and 
j , in the conservative case, let us now discuss the mod- 
ification of x brought by radiation-reaction. We define 
the supplementary contribution ^( RR ) to \ entailed by 
radiation-reaction by decomposing the total \ as 

X^ ot HE-J-) = X {consclv) (E-,J-)+X {RR) (E-,3-). 

(5.73) 

Here x ( ^ conscrv ' ) (-E'ii) i s function defined above in 
the conservative case and we have denoted by _E_ and j_ 
the energy and the angular momentum of the incoming 
state (considered in the infinite past, t — > —00). We 
are going to prove the following simple result concerning 
t^( rr ). When working linearly in the radiation-reaction 
J-'i, i.e. modulo terms that are formally quadratic in Ti, 
we can write 

x (RR) = ~ (x (conscrv) (^+^+) -x (conscrv) (£-,j-)) 

1 / r)v( conseTV ) a (conscrv) 

= I ( S^E+ ^ S^3 

2 V dE dj 

(5.74) 
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where S^E and S^j are the integrated losses of en- 
ergy and angular momentum, radiated (between t = — oo 
and t — +00) at infinity in the form of the corresponding 
fluxes and <&j. Note that (still modulo terms 0(J" 2 )) 
the result (|5.74p means that the total scattering angle 
^( tot ), in presence of radiation-reaction, can be written 
as 



X 



(tot) 



(E-,3 



■>-§( 



(coB S erv) ( £ +J+) 

X (conaeIv) (E-,j-)) . (5.75) 



Moreover, it can also be written (modulo 0(J- 2 )) as 

X (tot) (^-,J-) = X (cons ° rv) (^o, Jo) , (5.76) 



where 



E 



5* 

1 



E-) 



3o = -(J++3-) 



(5.77) 



are the average values of E and j over the incoming and 
outgoing states. As the radiation-reaction is of PN order 
T = 0(l/c 5 ), the accuracy of the results stated above is 
modulo corrections of PN order 0(l/c 10 ). 

To give a proof of the above statements, one should use 
the generalized method of variation of constants used in 
Refs. 26, 56, 61] , which considers the perturbation of the 
2PN accurate conservative dynamics by the radiation- 
reaction force. Moreover, one should extend the treat- 
ment of these references from the elliptic-like case they 
consider, to the hyperbolic-like one we are interested in 
here. This can be done, and yields a straightforward 
proof of the relations above. Here, for the benefits of 
simplicity, we shall content ourselves with presenting the 
proof of these relations in a simplified case where the un- 
perturbed dynamics is treated as being Newtonian, while 
the perturbing force Ti is considered at the fractional 
2PN accuracy. We shall, however, indicate the essential 
reason why the result still holds in the case where both 
the conservative dynamics and the radiation-reaction are 
treated more exactly, i.e. with a Hamiltonian of the type 

^(conscrv) = ^(Newt) + 1 ^(IPN) + }_ ^(2PN) > (5 7g) 



and a radiation-reaction of the type 



1 



T = J-(Ncwt) + jr 



1PN) 



J_jr(2PN) 



(5.79) 



When considering the simple case where the unperturbed 
dynamics is Newtonian, we can simplify the calculations 
of ^( RR ) by making use of the famous Laplace(-Lagrange- 
Runge-Lenz) conserved vector. Using scaled variables, 
f = r/(GM), j — J/(GMfi), p = p/[i (and, henceforth, 
dropping both the carets and the tilde's for easing the 
notation) we have the Laplace vector 



where j = r x p and n — r/r. Its time derivative is 
proportional to the perturbing force T (henceforth we 
shall also drop the tilde on F) and is given by 



dA 
~dt 



7xj+px (rx J) 



(5.81) 



If we write J~ in vectorial form, it has the structure 



T = a(r,p)p r n + /3(r,p)p 



(5.82) 



where the crucial information is that the coefficients a 
and j3 (which should not be confused with the quanti- 
ties introduced in the previous subsection) are time- even 
scalars, i.e., combinations of our usual scalars p 2 , p 2 , and 
1/r. [This holds for the 2PN-accurate a's and /3's.] In- 
serting this structure in the time derivative of A yields 

dA 

— = ap r ii x j + 2/3p x j . (5.83) 
dt 



Let us now decompose all vectors with respect to an or- 
thonormal basis e Xl e y , e z , with the x direction along the 
apsidal line (i.e. with e x a unit vector directed from the 
origin towards the periastron) and with e z being along 
the angular momentum: j = je z . We have 



n = cos (pe x + sin (pe y 

p = — [— sin 4>e x + (cos <f> + e)e y ] , 



(5.84) 



so that the two components of A = A x e x + A y e y read 
A x = ae sin 2 4> + 2/3(cos <j) + e) 



sin 4>{ae cos <fi + 2/3) , 



where we used the fact that 



p r = n • p = - sin ( 

j 



(5.85) 
(5.86) 



The crucial fact we wish to stress is that A x is an even 
function of <fi, while A y is an odd function of <j>. [Recall 
that the scalars a and j3 are functions of p 2 , p 2 and 1/r 
and are therefore even functions of <f>.] Remember that we 
have chosen the origin <po of </> at <po = 0, so that these 
parity properties of the vector A correspond to simple 
symmetry properties between the first half of the motion 
(between infinity and the periastron) and the second half 
(from the periastron back to infinity). When integrat- 
ing over time to get (at order 0{F)) the total radiation- 
reaction-induced change of A between —oo and +00, we 
deduce (using the fact that 4> — j/f 2 = (1 + e cos (j)) 2 /j 3 
is even in d>) 



6 (nn) A = A A 



(5.87) 



will be directed along the x axis. As the unperturbed A 
vector is simply 



A(t) = p x j - n (5.80) 



(5.88) 
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we conclude that the effect of radiation-reaction on A 
amounts to changing only the magnitude of the eccentric- 
ity e, without introducing any further angular rotation 
in the apsidal line. More precisely, as the magnitude of 
the perturbed A 2 (£) is given (at any moment) by 

A 2 (i) = P 2 j 2 + 1 - -j 2 = 1 + 2E(t)f(t) , (5.89) 
r 

where E(t) and j(t) are the instantaneous (Newtonian) 
values of the energy and angular momentum along the 
perturbed motion, we conclude that an incoming A vec- 
tor at t = — oo of the form 

A(t = -oo) = A = \Jl + 2E_j1e x (5.90) 
will end up, at t = +oo with the value 

A(t = +oo) = A+ = ^l + 2E +] \e x . (5.91) 

Let us now use these asymptotic results to compute the 
value of the scattering angle x^ tot \ including the cumu- 
lated effect of radiation-reaction. This is done by consid- 
ering the limits t — > ±oo in the defining expression (|5.80[) 
of A(t). Asymptotically, we have 



p(f = ±oo) = p± = ± J 2E + n± . (5.92) 



which is the relation that we have indicated above. 

Let us briefly indicate why this result extends to the 
case where the unperturbed, conservative dynamics is 
treated, say, at the 2PN accuracy. In that case one can- 
not use the Laplace vector because of periastron preces- 
sion. Instead one can use the version of the method of 



Here, C\(t) and 02(f) denote E(t) and j(t), while the third 
quantity c\ (t) corresponds to a possible additional angu- 
lar displacement of the apsidal line, beyond the effect 
linked to the radiation-reaction-driven adiabatic varia- 
tions of E(t) and j(t). The quantity c\(t) corresponds 
in our above simplified treatment to the direction of the 
vector A(t). We found above that the direction of A(t) 



Let us replace any vector V = V x e x + V y e y in the 
orbital plane by the corresponding complex number V — 
V x + iV y . In particular, the unit vector n(t) becomes 
the complex number n(t) = e 1 ^^. Its limiting values are 
n± = e lc ^ ± , where 0+ = 4>(t —> +00) and </>_ = <p(t —> 
—00). It is then easy to find that the asymptotic values 
of A(t) = A x (t) + iA y (t) are given by 

A_ = (-1 + i\j2E-j 2 _)n- 
A+ = -(l + iy/2E + jl)n+. (5.93) 
If we then define \± (and e±) by 

tan M = I = 1 (5.94) 
2 ^2E ±] \ ^el 1 

we conclude that 

A+ = -ie+e-^TH = -ie+e'^ e tlp + . (5.95) 

Our previous result show that A + has the same argument 
as A_. Therefore 

1 + ^ + *- = -%-% + *+ (5-96) 

so that the total scattering angle ^( tot ) = <j) + — <jy_ — tt 
(including radiation-reaction) is simply given by 

(5.97) 

I 

variation of constants used in Refs. [Ill [H^, [6l|, and 
adapt it to the hyperbolic case. Then the crucial quanti- 
ties which encode the effect of radiation-reaction on the 
scattering angle are the "varying constants" c\(t), 02(f) 
and c\(t) that enter the expression for <p(t) given in Eqs. 
(32b) and (33b) of Ref. (gy, namely 



(5.98) 

I 

did not include a secular change under the influence of IF, 
because of symmetry reasons linked, finally, to the time- 
odd character of J r . This fact has a correspondant in 
c\(t). Indeed, Ref. [6l| found that there were no secular 
changes in c\(t) (and ce(t)) precisely because dc\(t)/dt 
is an odd function of <fi, around the periastron, and re- 
marked that this was linked to the time-odd character of 



I 

x (consorv) (^-,.?-) + x (conscrv) (^+J+) 



I 

c/)(t) = f dt[l + k(c 1 (t),C2(t))]n(c 1 (t), 02(f)) + c x (t) + W(£; ci(t), c 2 (t)} . 

Jto 



24 



!F. When applying this result to a scattering situation, 
one again finds that the total scattering angle will be 
given by the average of the conservative ^-( consorv )(iJ , j) 
over the incoming (E-, j-) and outgoing (E +1 j + ) values 



of the two secularly-evolving "constants," E(t) and j(t) 
(i.e., Ci(t) and c 2 (i) in the notation of |61j|). 

Let us finally give an explicit estimate of the modifica- 
tion 



(5.99) 



of the scattering angle entailed by radiation- reaction. We 
will do this calculation at the leading PN order in the 
value of T, i.e., at the 0(1/ c 5 ) order only. We there- 
fore need the values of the losses of energy and angu- 
lar momentum during a hyperbolic encounter. From the 
(Newtonian-order) energy flux at infinity 



r 



we compute the integral 



dt$ E (t) 



(5.101) 



along the unperturbed motion, using cf>, rather than t, as 
integration variable, i.e., 



5Gc 5 \ r J 



11 



S^ RR) E = E + -E- 



(5.100) Computing this integral, we find 
I 



GMfij 



$ E - (5.102) 



2v 



15c 5 j 



5„-7 



-(673e^ +602)Je 2 _ 

o * 



1 + (37ei + 292e 2 _ + 96)^ (e_) 



(5.103) 



where </>5j_(e_) is defined (in keeping with Eq. (|5.53[) ) as the Newtonian angular momentum flux at infinity, 



6^_(e_) = arccos ( j = — + arcsin ( — 



(5.104) $ J= ^ {jG M)(^Y (2v 2 -3r 2 + ?\ (5.105) 



This result agrees with Eq. (2.10) in [62| . Similarly, from we computed 
I 



§ (RR) 



3=3+- 3- 



0+ 



0+ (GM) 2 vj M 



(5.106) 



We find 



5 m) ] = 



(2ei + 13) Wei - 1 



5c 5 ji 
+(7 e 2 _+8)^(e_)] , 



(5.107) 



I 

of E and j, we are mainly interested in the radiation- 
reaction-driven change in the eccentricity, namely 

sm , e = de(E,j) giRR)E de(E,j) s(KR) 
dE dj 



where <^(e_) is the same function as above. Using the results above for S {R ^E and <5( RR )j we find 

As the (conservative) scattering angle is a function of 2 e 

the eccentricity, i.e., the combination 5' RR ^e = ^-Q(e-) , (5.110) 

15^ c°j_ 



e(E,j) = Vl + 2£j 



(5.108) where 
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Q(e-) 



- 1 



3e 2 



■(72ei + 1069c 2 . + 134) + (304 + 121e 2 _)<f>%(e-) . 



(5.111) 



We have also checked this result by computing the change 
in the Laplace vector A. We find that the </>-derivative 
of the associated complex quantity A = A x + iA y reads 

^=*vje i + [-3yV 8 + 6jW a 
aq> 15c D 

+i(7f - I5j 2 )r 2 r' - 12(f + j 2 )f 3 ] , (5.112) 

where the prime denotes a (^-derivative. Inserting the 
Newtonian orbit f = j 2 /(l + e cos cf>), and integrating 
between 0_ and <j)+ yields 

S m) A = A+ _ A _ = _2 -|=g-Q( e _), (5.113) 
15 c J jl 

in agreement with Eq. (|5.110[) Finally, as 



X (conscrv) (e) = 2arccos [ - 2 j - tt 



= 2arcsiu | - ) (5.114) 



we have dx {consci ' v) / de = -2/(e\/e 2 - 1) so that 



2 ae 
j(RR) e 



(5.115) 



Finally, the radiation-reaction contribution to the scat- 
tering angle is given by 



where Q(e_) is defined in Eq. (|5.111|) . 



VI. SUMMARY AND OUTLOOK 

Let us summarize the main results of our work: 

1. We have introduced a new approach to the compu- 
tation of the gravitational radiation-reaction, based 



on the identities (|2.15|) , (|2.16j) satisfied by the com- 
bined energy and angular momentum flux function 
Eq. (233- 

2. We have computed some "minimal" version of 
the 2PN accurate radiation-reaction force .F(x, p) 
which must be added on the rhs of the Hamil- 
tonian EOB equations of motion when describing 
general orbits (elliptic-like or hyperbolic-like) . The 
radial, J- r , and azimuthal, J-^, components of the 
radiation-reaction force are explicitly given as func- 
tions of the EOB position and momenta by Eqs. 
(|3 . 62[) and fl3.63p . Our calculations were based 
on the transformation properties of the three ba- 
sic scalars X\ ~ p 2 //i 2 ~ v 2 , X 2 ~ P 2 /^ 2 ~ r 2 
and X3 ~ GM/r between the various coordinate 
systems used in PN theory (harmonic, ADM and 
EOB). 

3. We have also computed the "Schott" contribution 
to the energy, corresponding to the above minimal 
construction of J- . It is given as a function of the 
EOB position and momenta by Eq. f|3 .69[) . In par- 



ticular, we pointed out that E, 



(schott) 



does not van- 



ish during quasi-inspiral but is proportional to p r 
and is given by Eq. (|5.7|) . 

4. We provided a new understanding of the gauge free- 
dom in the construction of the radiation-reaction. 
It is linked to the arbitrary choice of (i) the Schott 
contribution to the angular momentum, and (ii) the 
part of the Schott energy which is proportional to 
the cube of the radial momentum p r . This explains 
very simply why there exist 2x1 arbitrary param- 
eters in T at the Newtonian order, 2 x 3 at the 
1PN order and 2 x 6 at the 2PN order [and then 
(n + l)(n + 2) at n PN order]. 

5. We pointed out that there is an inconsistency be- 
tween the assumptions that are standardly used in 
current implementations of the radiation-reaction 
force in the EOB formalism, namely Eqs. (|5 . 8[) and 
(|5.9[) . We showed that if one adopts the assump- 
tion J-0 = —$7 (which is convenient, and always 
possible) this essentially determines (during inspi- 
ral) a nonzero value for the radial component of 
the radiation-reaction force, given by Eqs. (|5.14l) . 
dOffl and ([BTTTp . 

6. We introduced a new way of parametrizing (con- 
servative) hyperbolic orbits in PN theory, by the 
simple quasi-conic equation (at 2PN) 
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E + i_ a2 C P) 2 + _L a3 C P") 3 + _L a4 C P") 4 = i 

r c 2 V r / c 4 \rJ c 4 V r / 



A" 



(6.1) 



and emphasized that the two quantities e ("eccen- 
tricity") and K ("periastron advance") are gauge 
invariant. The gauge-invariant eccentricity e is re- 
lated to the scattering angle x an d to K via Eq. 
(|5.40l) . Moreover, K and \ are given, in EOB the- 
ory, by simple (complete or incomplete) integrals 
over the inverse-radius u = GM/r, Eqs. ()5.64l) . 

7. We have showed how the effect of radiation-reaction 
on the scattering angle can be computed (modulo 
correction 0{T ) = 0(l/c 10 )) from the sole knowl- 
edge of the losses of energy and angular momen- 
tum at infinity, see Eqs. (|5774"|i . (|5775|l and ([5775)1 . 
This result might be used to subtract the effect of 
radiation-reaction on the scattering angle obtained 
in numerical simulations, by using only numerical 
data in the asymptotic domain at infinity. We also 
gave an explicit expression, at leading order in 1 /c, 
for the additional contribution to the scattering an- 
gle due to radiation- reaction, see Eq. (|5.116[) . 

Finally, let us point out some of the future research 
directions that would complete our results: 

(a) In the present work we have not included the effects 
of tails on the radiation-reaction. We plan to treat 
this issue in a future publication. 

(b) Here we obtained the components of the radiation- 
reaction force !F in the form of a standard, non- 
resummed PN expansion. However, the cur- 
rent most successful implementations of the EOB 
formalism make a crucial use of efficient re- 
summations of in the circular limit. It would 
be interesting to concoct resummation schemes in 
the more general context considered here. For in- 
stance, in the case of slightly elliptic orbits one 
might hope to improve the numerical validity of 
our PN-expanded FiS by first factorizing the "cir- 
cular part" of these components, and re-summing 

I 



them by the method introduced in [14|. We gave 
some partial results towards this goal in Sec. IV. 

(c) Let us finally mention that, in order to have a com- 
plete EOB formalism for general orbits, there re- 
mains the problem of expressing the emitted grav- 
itational waveforms in terms of the EOB phase- 
space variables. The transformation formulas we 
provided should be also useful in this respect. 
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Appendix A: The 2PN energy and angular 
momentum fluxes in the far zone in EOB 
coordinates 



1. Energy flux 

The 2PN energy flux (excluding tail terms), scaled as 
in Eq. (|3. 281) . can be written as 



&* oh \x\) = [XIY {C A X% + e 2 C AB X%X° B 



-e C ABcX\X e B X e ^) 



where 



C A = -v 
5 



(Al) 



(A2) 



and 



3^^(1347- 2556^) 



Cii 

C 22 = -o|f^(-18549 + 108i/) 



C12 - ^^(-1333 



C 2 



9^(12960 



-412i/) 
2304z/) 



Gig = -^K36720 ~ 6696^) 
C 33 = 2^(432 - 1728*/) 



Cm 

Cn3 
C.23 
C222, 
C233 



= ^i/(-159 - 838i/ + 1874i/ 2 ) 
= -^v{12,9Qv 2 - 6362i/ + 4761) 

= _4_;y(_99/19,/ _|_ 81 1 1 _L IRR,,^ 

= 315K- 

2 v(7MAv 2 - 32052j/ + 81263) 



2242z/ 
4498t/ 



8111 + 368i/ 2 ) 
2828i/ 2 - 2501) 



2835 



C112 
C122 
C133 
C223 
C333 



3^i/(490t/ 2 + 1523i/ + 1101) 
-j|gi/(-372i/ + 2034i/ 2 + 43) 
gjjgz/(22836i/ 2 - 61596i/ + 198961) 
2^5i/(590i/ - 58611 + 974z/ 2 ) 
-g^z/(145- 648i/ + 272i/ 2 ). 



(A3) 



(A4) 



r 



Note that we are listing here and below the indepen- 
dent components of the symmetric "tensor" Ca 1 ...a„- 



When explicitly effecting the multisummations present 
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in the contractions Ca 1 ...a„Xa 1 ...a„ (with Xa 1 ...a„ = -^Ai • • • -^A„) they appear multiplied by the symmetry 

factors of Eq. p. 251) , namely 



CabXab = CuXti + C21X12 + C33X33 + 2C\iX\2 + 2C13A13 + 2C23A23 
CabcXabc = Ciii-Xiii + C222X222 + C333A333 + 3C112X112 + 3C113X113 

+3Ci22^122 + 3C133X133 + 36*223^223 + 3C*233^233 + 6C123A123 . (A5) 



I 

2. Angular momentum flux where 

The 2PN angular momentum flux (excluding tail 
terms), scaled as in Eq. (|3. 281) . can be written as = —v[2 —3 2] (A7) 



5 



^j° h) (X A ) = J [X!] 3 (B A X A + e 2 B AB X%X% 

+e i B A BcX A X e B X e c ) (A6) 



B u = 3^^(330 - 1272j/) B l2 = ^^(-396 + 900;/) B 13 = ^^(-5928 - 3288?/) 

B22 = -3I5K- 1710- 1080;/) B 23 =-^^(-11520-6001/) £33 =- ^(11898 - 1548;/) 



(A8) 



and 



S m = 3^i/(-1051z/ - 50 + 750i/ 2 ) B112 = ^sK 205 ^ 2 + 3347?/ - 971) 

fl U3 = g|g^(7802z/ 2 - 6057 + 6238?/) B122 = -^;/(872;/ 2 - 430 + 1489i>) 

B123 = -jgw i/(5516i/ 2 - 10392;/- 26869) £133 = -Jgi/(4312 + 4448i/ 2 - 21843;/) (A9) 

-B222 = |j/(76z/ 2 - 78 + 155*0 B 22 3 = ^v{272v 2 - 2725;/ - 5255) 

B233 = gfgi/(-2997 + 1098;/ 2 - 6566;/) £333 = ^^(46085 + 3042z/ 2 - 6741;/) . 



3. Combined energy and angular momentum flux 

The combined flux $£° b) = $£ ob) - GM0 e $^ oob) (ex- 
eluding tail terms) can be written as 

^t\x%) = [xi\ 3 Q 2A {x%) (A10) 



where 



and 



Qn 
Q22 



24 



Q 



!12 

23 



= Av 



Ql3 = p 

Q 33 = o 



(All) 



Qui =-^(-174- 17761/) 
<3ii3 = 45^(10134-2328z/) 
Q123 = g|5J/(-18234 + 1536i/) 
Q222 = -5TBf(-1710 - 1080;/) 
Q233 = -3M5^(-12510 + 1548i/) 



Qii2 = -45^(534 + 3432^) 
Qi22 = -gi5^(534 + 3432;/) 

Q133 = 9I5-K-3690 + 468;/) (A12) 
Q223 = -j^s ^(-76194 + 2952;/) 
Q333 = j&g"^ - 1728;/) ; 
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finally, the 15 independent components of Qabcd are 



Q1111 = -3^5^(163 - 202v + 1764i/ 2 ) 
Quia = I i 5 i/(456i/ 2 - 16682i/ + 1905) 
Qnas = i§UKl4159z/ - 5543 + 7784^ 2 ) 
Q1222 = -^v(6916v + 3788z/ 2 - 2211) 
Q1133 = g|5^(790i/ 2 + 3200z/ - 19679) 
Q2223 = 2io^(-6779i/ - 9211 + 2428j/ 2 ) 
Q3333 = 1^(1768^ - 18361/ + 24455) 
Q3332 = -5^o^(5400^ 2 - 40068;/ + 118193) . 

Similarly to Eqs. (|A3I) and (|A4[) above, the symmetry factors multiplying the independent components of the sym- 
metric tensor Qabcd are given by 

QaBCdXabCD = Qllll^llll + <52222^2222 + <53333-^3333 

+4<5lH2^1112 + 4<5lll3^1113 + 4(51222-^1222 
+4Ql333^1333 + 4Qa223^2223 + 4(52333-^2333 
6Qll22^1122 + 6(51133X1133 + 6(52233-^2233 

+ 12(5n23A"ii23 + 12<5i223^1223 + 12(5l233^1233 (A-14) 



Quia = -^go^(3372i/ 2 + 9424*/ - 3445) 
Q1122 = ^1/(10984*/ + 6252j/ 2 - 2959) 
Q1233 = -1^0^(3296^ - 2134*/ - 108839) 
Q1223 = - 3^o^(-30434z/ - 66503 + 30312i/ 2 ) 
Q2222 = §K76^ 2 - 78 + 155*/) 
Q2233 = 5M5^(3991^ 2 - 9005*/ - 133530) 
Q3331 = t!oK1768z/ 2 - 1836*/ + 24455) 



Appendix B: Hamilton equations in EOB 
coordinates: expansion at 2PN 

• Equation for r e 

From Hamilton's equations we have 



dH(cob) 

dpl e) 



with 



where 



C(n e ■ p e ) = Cp^ (Bl) and 



C = 1 + ^CaX'x + £ A CABX e A X e B , 
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(B4) 



• Equation for pi 

From Hamilton's equations we have 



In X\ variables we have 



an 



(eob) 



dr K 



%^ = xic lt3 (x%) 
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and 



(B5) and 
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where 



CU = [1,-1,-1] 



(B7) 



Appendix C: Schott energy at 2PN in EOB 
coordinates: minimal gauge expression 
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The minimal gauge expression for the Schott energy is given by the first of Eqs. (|3.69p . that is 



E. 



schott 
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Appendix D: Radiation reaction force at 2PN in 
EOB coordinates: Minimal gauge expressions 



where 



The minimal gauge expression for the radial compo- 
nent of the radiation reaction force is given by the second 
of Eqs. (|3T69|) . that is 
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The minimal gauge expression for the azimuthal com- 
ponent of the radiation reaction force is given by the third 
of Eqs. IglgD , that is 



r. 



(cob) _ J r v e]3 



±{xiys h3 (x_ 
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and 
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Finally, for the expression Eq. (|3.70p of J>(X|,X|,X|) in terms of the EOB variables X|,X|,X| we have the 
following coefficients 
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Similarly, for the expression Eq. Q3.70P of T^iX^, Af, X|) in terms of the EOB variables (having also 

used the expression (|3.7ip for j), we have the following coefficients 
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Appendix E: Coordinate transformations in phase 
space: harmonic, ADM, EOB 

In PN theory there exist at least three different co- 
ordinate systems that are largely used: harmonic (h), 
ADM (a) and EOB (e). Each of these systems has its 
own utility and we shall discuss here their transforma- 
tion laws at the 2PN order. We work with the scaled 
position variables q/j = x/,/(GM), q a = x a /(GM), 
q e = x e /(GM) and similarly for velocity or momentum 
(per unit reduced mass) variables, which are simply de- 
noted by p/i, p a , p e without recalling the tilde notation. 



Phase space variables associated with harmonic coordi- 
nates are only (q^, v^) (no ordinary Hamiltonian exits in 
this case), whereas for the ADM and EOB cases one has 
either (q Q ,v Q ) and (q e ,v e ), respectively or (q Q ,p a ) and 
(q e , p e ). With each choice of phase space variables (h,a, 
or e) is associated a family of fundamental scalars, that 
is for example 

(qfc, vfc) -4 X[ l = v 2 , A 2 " = (n h ■ v,,) 2 , X% = — 

where n h = q.h/<lh> e ^ c - We list below the main transfor- 
mation laws among phase space vectors as well as funda- 



mental scalars. 

J 



1. ADM vs harmonic coordinates 

ADM vs harmonic phase space vector 2PN-transformations are the following: 

1) (Qo.Vo) -> (fVi,\h) 
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Concerning the transformation of fundamental scalar quantities, we recall once more the notation introduced in 
Sec. II, namely X\ — v\, X 2 = (n^ • v^) 2 , X 3 = The same notation for the ADM variables leads, as explained 
before, to the two possible choices 

1 

la 



X?=p 2 a , X 2 Q = (n a -p a ) 5 



and 

We find explicitly 
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Equivalently, using our "tensorial" notation Eqs. ()E9[) are summarized by 
1. y« = iQ?^^), with iQl' 1 = 5ib and the iQ^ all vanishing, while 
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Similarly, we may summarize the relations X\ = ^Cf 3 (X^) as indicated below 
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2. X 2 ° = 2 C$(X%), with 2 C« h - 6 2B , 
and finally 
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3. X 3 Q = 3 Cf h 3 (XX), with 3 C^ 1 = S 3B ; here all the 3 Cf B vanish while 
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2. Harmonic vs EOB coordinates 

Harmonic vs EOB phase space vector 2PN-transformations are the following: 

1) p e «• V e 



v e = q e = p e + e 
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We recall that p e denote momenta per unit reduced mass (indicated without the tilde, for convenience). 
2) (q/nVfc) «-> (q e ,p e ) 
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+e 4 



i/(lli/-3) 4 ^(19 + 9^) 2 3z/(5z/-3) 

o V h V h - 



8q h 



, 2 1/(1/ -19) 



q^ 



'v(7v-l) 2 3j/(3 + 5j/) . 



Vh + e 



3^ + 2 
2g h 



4% 
(n,j ■ v h )n h 



\-2v 3^(3 + 5i/) 

vi + 1 v h 



+e 4 \ (n h ■ v h )n h 



(2lv{l + v) 



(n/j • w h ) 2 - 



23^ 2 -15^-4 9 15i/ 2 + 23^+12 



-vt + ■ 



+ 



- - 2i/ + 3^ 2 w,t 



21i/ 2 + 7^-24 2 v{l7v + 25) 



+ - 



v 2 - 14i/ + 3 



2^ 



(E20) 



q^ = q e + e 



:P~ 



i/ + 2 



2 e 2<7e 

i/(5i/ + 17) 



q e + i/(q e • p e )p e 



+i/(q e • Pe 



8<Ze 
1/ - 1 



2 1/(1 + I/) 4 , l/(3l/ - 1) 2 1/(1/ -19) 
(n e • Pe) Ve + — P, 



Vh = Pe 



+ £* 



2 

3^ + 2 



-Pe + 



v - 19 



4<? e 

(n e • p e )n e - 



v + A 2i/-l 



+ 



"Pe Pe 



2<z e ^ — c V 2<z e 2 
(Zvlbv + l), N2 7i/ 2 + 23z/-4 2 3^ 2 -9iy-4\. 

7, ( n e ' Pe 7. Pe "j 9 1 I 11 * ' Pp ' n '- 

V 8g e 8g e Aqj 
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,. 15i/ 2 - 29i/ - 8 , x2 3 -81/ 4 7v 2 -41v + 8 2 v 2 -\hv-\ 
+ ( 8^ (ne ' Pe) + + 8^ Pe 2^ 



(E21) 



3. X§ = 3 Cf%(X>A), with 3 Cf = 5sb 

rie. - 
3^13 



and finally 



Ceh 3,, 2 i v rieh v i 1 ,, rieh 5 ,.2 , 11,, 

• 113 - + 8 3 °123 - ~— + V2 V 3L/133 - + 

3C223 = T" 3CI33 = f!^ 2 " 2^ 3C|33 = ^T v + 1 ■ 



Let us consider now the transformation law of the fundamental scalars X\ and X\ as represented by X\ = 
A C^ 3 (X%), with 

1. = iC^kxd, with = 

i<7fr = -l + 2i/ iC& e = -|-2 iC^^-fiz-l, (E22) 

and finally 

rihe ,,2 q,. 1 1 /"i/ie v 55.. 1 4 rihe 17 ,,2 1 1 ,, 1 2 

1°111 - v + 1 iO n3 - — - j^V + 3 1^123-^12^ +6^+3 (F?^ 

r»/ie _ !/ 2 1 17,., 5 r , ' le —5,. 2,1,. rihe _3,,2 , 29,. , 7 V^^-V 

1^133 - + T ^+ 3 lG 223 - i V + 4 V 1^233 — ~i v + ~ v + 3 ■ 

2. ^ = 2 Cl%{X%), With 2 C* e = 6 2 B, 

2 C[ h % = -\+2v 2 C% = -2v 2 C^ = -2>-2v, (E24) 

and finally 

2 C^ 2 = 2v 2 -2v+\ 2 C^ 3 = -^-fu + l 2 C%% = l^ + ^+l 

2 Cf| 2 = -3^ 2 + z/ 2 C% 2 =4v 2 2 c% % = ^ + fu + 4. ^> 

3. X$ = 3 C^ 3 (X%), with 3 C^ = d 3B , 

3^ = -! 3^1 = -! 3^1 = 1 + 1, (E26) 

and finally 

rihe v 2 < v rihe u rihe 7 , ,,\ 

3^113 — — + 24 3^123 — 12 3^133 — "SI^ 1 + fE27l 

rihe _ v rihe 5 ,.2 , 13,. rihe _u 2 15,,,-\ ^ ' 

3^223 - T 3^233 — "24^ + "§" 3^333 - T _ + 1 • 

Similarly, for the transformation X\ = aC{ h 3 (X\) we have 

1. X\ - iCf%(X'X), with x C B h = S 1B , 

!Cft = l-2v 1 C 1 t = | + 2 iC|J = |i/ + l ) (E28) 

and finally 

iCft = 7v 2 -bu + l X C\\ 3 = -f 2 v 2 iC$ 3 = -fv 2 - \v + i 

i^f& = ^-4i/+i 1^ = ^(1 + ^) iC 2 t 3 = f- 2 + f-+f • ( ^ 9) 

2. Xf - 2 Cf h 3 {X%), with 2 <7f, fc = S 2B , 

2C12 =5- 2 ^ 2C 2 t - 2^ 2 Cf 3 l = 3 + 2^, (E30) 

and finally 

rieh _ c.,2 8„, 1 ne* — 53 ,,2 29,. ,5 ,oefc —25,2 ,23,, 

2Wl2 - bly ~3^+3 2^123 " -T2" _ T2^+3 2^323 " + — ^ /poi\ 

2<7& = -51/2 + 1/ 2CI 2 l 2 =4^ 2 ' 2C|& = f^ + $i/ + 6. 



°13 - 4 3<^ 23 - 2 3^33 - ~ 2 _ 1 ' 



(E33) 
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3. EOB vs ADM coordinates 

EOB vs ADM phase space vector 2PN-transformations (see e.g., sec. VI, Eqs. (6.22) of ref. [![) are the following: 
1) (qo,v ) O (q e ,p e ) 



q e = q a + e 



+ - 5-- + - 



2/ g a 



+ (Qa • Pa) 





1 " 




(>-5) 




Pa 




9a _ 





Pa 



(^•P a )^(l+2)^+2^-^( 1+ 2) Pa 



£(io - „)rf + Us + + h-2 18, + , 2 )1 

8 8 q 4 g a 

1 



+e <{ (q a • p a ) — 



(E34) 



q a = q e 



^^z( 1+ D) qe+Kqe - pe)pe 



^{[-5(1 + ^ + 1^5-1 

+(q e • Pe 



A&_ v ( 2+ *„) (qe-Pe) 2 + ^ 
9e V 8 / of 



7,-1 1 



,(, — 1 o VI v\ 1 



2 ^ ' 2 



9e 



Pe 



Pa = Pe 



(1+ 2J(qe-p e )^q e + 



-e 4 i (q e 'Pe)-, 
' 9e 



3 (v ,\ 2 ^ 3 2 (q e - Pe ) 2 / 3 5 



3 A i 



2 2 A J q e \ 

4 2 J g2_ 



(E35) 



the transformation law of the fundamental scalars X\ and X\ as represented by X\ = aTi%{X\ 



Let us consider now 
with 

1. X? = X T^{X%), with x Tg e = <5 1B , 

iTff = -v !Tf 3 e = 1 + | iT; 
and finally 



'ae 
23 



Tae _ ..2 j_ 1 ,, 
. Ill - ^ + 4^ 

, rpae 5 , .2 1 _i_ 7 

lJ 133 ~ v2 V _ 6^ + 6 



l7iT 3 = -P-2^ 1^3 = ^(1 + ^) 

7^ae iv 2 rpae 7 T ,2 i 4 T , 

1 J 223 — X 1 J 233 — If" ~r 3 



2. X 2 a = 27 1 r 3 (Vl), With 2 T° e - ,525, 

and finally 



pae v rpae o,/ 

2-i 12 — 2 2J 22 — 



2 J 112 

T'ae 
_ 9' 



2^122 

T'ae 
- 222 



-V 2 + \v 

4, 2 



T>ae 1..2 1 . . 

2^123 — 1 V ~ 2 V 
~<ae 1 ,,2 i 7 



2^3 = - -'' 



(E36) 



(E37) 



(E38) 



(E39) 
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3. XI = 3 T%{X%), with 3 T» e = 5 3B , 

^ = -7 ^§ = -^ 3 T£f = | + l, (E40) 



and finally 



rpae v i v rpae v rpae 7 .,2 f 

3 J 113 — 8j 24 3^123 — 12 3 J 133 — 2 24 " 12 

T>ae f o^ae 5 ..2 _i_ 5.. rpae v 3,. _i_ 5 

3^223 -— 3^233 ~~2i V +6^ 3 ^ 333 - — ~ 4 ^ + 4 



and finally 



rrea ,,2 f ZTea 5 , .2 11,, rrea f f 

l-^lll - V ~J 1^113-^12^ _ ~ V ltL 12Z - ~ V2 

Tjea v 2 I 13., i 5 Ilea 3 , ,2 i a,, XJea v 2 10,, 2 

1^*133 - 12 + + 6 1^223 - 4^ + lV 1^233 — 12 3~ — 3 



(E41) 



Similarly, for the transformation law = aH^^X^) we have 
1. X{ = !Hf%(X%), with ii/|f = 5 1B , 

iH{l = v x Hft = -\-\ i^2 e 3° = f + l 5 (E42) 



(E43) 



2. Xf = 2H?%(X%), with 2 H% e = 5 2B , 



and finally 



Tjea 
2^12 



'22 



2v . 



(E44) 



rrea v 

2-^112 — 4 

rrea 12 5. 

2-^223 — 6 V ~ 3 l 



Tjea ,,2 

2-"l22 — ~ v 

2 HII 2 = 4^ 2 



-^ 2 + 



'123 — 12* 

rrea 1 ,,2 7 . 1 1 

2^233 -fi^ -6^+6 



(E45) 



3. Xf = 3 tfr 3 (*2), with 3 T» e - 5 3B , 



Tjea 
3^13 



rrea v 

3 H 23 - 2 



rrea v \ 

3^33 - _ 2 ~~ 1 • 



(E46) 



and finally 



rrea _ 1/ 
3^113 - "8" 

rrea _ v 
3-^223 — T 



24 



rrea 

3^123 - 

rrea 

3-^233 — 



3 



rrea _ v 
3^133 8" 



rrea 11 ,, 1 

3-^333 - T + 



(E47) 



4. Transformation of the angular momentum 
variables 



While the conserved angular momentum of the system, 
J, has its usual, simple expression in ADM and EOB 
variables, namely (in reduced form j = J//xM) 



J = q a x p a 



q e X Pe . 



(E48) 



its expression in harmonic variables involves an extra PN- 
correcting factor f h = l + (3(l/c 2 ), namely 



j = fh Qh x v h . 



(E49) 



where 



f h = l + e 2 (-1(3* - 1)X? + {v + 3)X^j 



-(4^ + 14-41^ 3 



-l(9z, 2 + Wu - 7)X{ 1 3 - 1(5^ 2 + 2v)X* 3 



+ |(1 + 13^ - 7v)Xh 



(E50) 



Appendix F: Some reminders of Newtonian theory 

The relative motion of two bodies with masses mi and 
m 2 can be treated as that of a single body with effective 
mass fj, = mim 2 / (mi + m 2 ). Indeed, after separation of 
the motion of the center of mass (with M = mi+m 2 ) 



R = 



mixi + to 2 x 2 
M 



(Fl) 
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one gets the following Lagrangian for the dynamics of the 
relative motion 



,'l.o GM 



(F2) 



where r = xi — x 2 and r = |r|, from which follow the 
momenta 



p = fir = fix 
and then the Hamiltonian 

/ p 2 GM 



(F3) 



(F4) 



We systematically use a "tilde notation" for quantities 
per unit reduced mass; for example 

Co = Co/ ft , Ho = Ho/fi . (F5) 

The conservation of the angular momentum 

J = r x p = fir x v = /iJ , (F6) 

allows one to study the motion in the x — y orbital plane 
(orthogonal to J = Je z ). Using polar coordinates x % = 
(r, (j)) leads to the Lagrangian per unit reduced mass 

Co(r,r,^j>) = ±(f 2 +r 2 <j> 2 )+™, (F7) 



so that 



dC 



dCo 



or d(f> 



and 



GM 



(F8) 



(F9) 



The dynamics simplifies if we use the following rescaled 
variables 



J . t 
t 



- - r ~ - El ■ - p <t> 

r ~GM' Pr ~ fi' 3 ~ GM ~ GMfx ' " GM 

(F10) 

The Hamiltonian corresponding to these scaled variables 
is 

1 f~2 , f\ 1 



Ho(p r ,r,j,<f>) = -[p r + — 



(Fll) 

2 \ r z J r 
and the equations of motion read 

!tt-j_ dp\_f }_ , P19 x 

The integration of the radial equation fully determines 
the orbit 



f (</>) = 



1 + eo cos <f> ' 



P = ] 



(F13) 



or 



f(4>) 



= 1 + eo cos < 



also implying 



df e . 
— = — sini 
dt j 



eo being the eccentricity of the orbit given by 
e (E,j) = l + 2Ej 2 ; 



(F14) 



(F15) 



(F16) 



where E = Ho is the conserved energy per unit reduced 
mass. 

One has now to distinguish among the various types of 
orbits: elliptic (0 < e < 1; e = in the circular case), 
parabolic (e = 1) and hyperbolic (eo > 1). 

• Elliptic orbits The solution of the equations of mo- 
tion can be given in terms of the eccentric anomaly 
u as follows 



n(t - 1 ) 
f 



u — eo sin u , 
fio(l — eo cosu) . 



2arctan 



l + e 

1 - e 



tan 



(1) 



where 



}3 ' 



do 



1 

2E 



(F17) 



(F18) 



do being the scaled semimajor axis of the ellipse, 
«o = do /{GM). Other useful relations are 



J = 



implying 



y/-2E ' 
P = f = oo(l - eo) : 



(F19) 



(F20) 



and 



x = f cos(</> — 0o ) — a (cosu — e ) 
y = f sin(0 - <po) = d ^Jl - e\ sin^ 



(F21) 



The circular orbit case, i.e. r = = r, corresponds 
to eo = 0, that is 



f=f, E 



2f 



or 



1 + 2Ef = . 



(F22) 



(F23) 
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• Parabolic orbits The parabolic case (E = 0) is ob- 
tained from the elliptic one taking "consistently" 
the limit E 0. For instance, in Eq. (|F17j) one 

poses 



\J -2Ex 



(F24) 



and takes the limit E — > keeping x fixed. The 
result is 



t-tn 



2 
2.r 



(F25) 



• Hyperbolic orbits Transition to the hyperbolic case 

(F26) 



is accomplished by the substitution 



in the elliptic case relations, so that 

n(t — to) = — u + e o sinh(u) 

f = fio(l — eo cosh(u)) 



= 2arctan 



eo + 1 , , (u 

tanh — 

e - 1 V 2 



with 



1 

2E 



(F27) 



(F28) 



The "parameter" p entering the polar form of the 
orbits is still given by 



P = f = a (l ~ el) . 
The scattering angle is given by \54 



, X 
tan — = 

2 



1 



1 



l 



(F29) 



(F30) 



2Ef 



where eo = y 1 + 2Ej 2 . Note also the equivalent 
relations (whose 2PN analogs we often use in the 
main text) 



X ( 1 

— = arccos 

2 V e 



1 



arcsm | — 

,eo 



(F31) 



and p (max ) = P0/r min = jGM/r m[n = j/f {min) . 
Hence, 



J ^*(min)P(max) ) 
1-0 1 



E — 2^(max) 



(F32) 



so that 



1 + 2Ej 2 = (?5 2 max) f (min) - l) 2 , (F33) 



which can be replaced in Eq. (|F30[) if one wishes 
to express tanx/2 in terms of P 2 max) f( min ). 

Anticipating applications of our framework to numer- 
ical relativity simulations of hyperbolic encounters, let 
us indicate an estimate of the simulation time t( stop ) 
(counted from the periastron passage) necessary for ex- 
tracting from the corresponding polar angle 4>+(t (stop)) 
(counted from the periastron) the scattering angle x with 
some prescribed accuracy e = 10~ N <C 1. 

Consider the Newtonian relations for hyperbolic mo- 
tion with to = = <f>o, i.e., 



eo sinh u — u = nt 
tan 



eo + 1 



eo 



- tanh (|), (F34) 



where n — |&o| 3 ^ 2 . The asymptotic value for cf> corre- 
sponds to u — > +oo, that is 



tan 



e + l 
2 ) V e - 1 

From Eq. (|F35|) also follows 
tan <j) + 

so that 

tan (i: 



el-l: 



1 



tan( 



(F35) 



(F36) 



(F37) 



Let us define an "incompleted" instantaneous scattering 
angle x(t) by 



* , x(t) 



(F38) 



From Eq. (|F34|) . x(t) satisfies (when it is large and pos- 
itive) 



Xc 



cot ( — j = cot j I '.null — = cot (J^^j T 

Xo 



u 
2 



Xoo\ l-e-° 



cot(^)(l-2e-«) 



(F39) 



The scattering angle can also be expressed in terms 
°f ^(min) an d P(max)- Indeed, at the point of mini- 
mal distance (periastron) r — ?"( m i n ) one has p r = 



cot (|) 

V2/ ~l-2e~ u . 



cot (*f ) 



(F40) 
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From the "time" equation (IF34|) i evaluated for large u, 
i.e., 

e y«ni (F41) 

we have 



The condition for the left hand side of Eq. (|F43|) to differ 
from 1 only within some precision e = 10 — is then 

10- w « (F44) 

^*(stop) 

that is 



so that 



cot (2) 



cot I 



go 
2n£ 



1-^ 



nt 



(F42) 



(F43) 



c (stop) 



1 + 2Ef 

(2S) 3 / 2 



■10 



A' 



(F45) 
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